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Abstract 

We derive probabilistic limit theorems that reveal the intricate structure of the phase 
transitions in a mean-field version of the Blume-Emery-Griffiths model [4]. These prob- 
abilistic limit theorems consist of scaling limits for the total spin and moderate deviation 
principles (MDPs) for the total spin. The model under study is defined by a probability 
distribution that depends on the parameters n, (3, and K, which represent, respectively, the 
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number of spins, the inverse temperature, and the interaction strength. The intricate struc- 
ture of the phase transitions is revealed by the existence of 18 scaling limits and 18 MDPs 
for the total spin. These limit results are obtained as K) converges along appropriate 
sequences {f3n-,Kn) to points belonging to various subsets of the phase diagram, which 
include a curve of second-order points and a tricritical point. The forms of the limiting 
densities in the scaling limits and of the rate functions in the MDPs reflect the influence 
of one or more sets that lie in neighborhoods of the critical points and the tricritical point. 
Of all the scaling limits, the structure of those near the tricritical point is by far the most 
complex, exhibiting new types of critical behavior when observed in a limit-theorem phase 
diagram in the space of the two parameters that parametrize the scaling limits. 

American Mathematical Society 2000 Subject Classifications. Primary 60F10, 60F05, Sec- 
ondary 82B20 

Key words and phrases: scaling limit, moderate deviation principle, second-order phase transi- 
tion, first-order phase transition, tricritical point, Blume-Emery-Griffiths model, Blume-Capel 
model 

1 Introduction 

The purpose of this paper is to analyze a new set of phenomena associated with the critical 
behavior of probabilistic limit theorems for a mean-field version of an important lattice-spin 
model due to Blume, Emery, and Griffiths [4]. These probabilistic limit theorems consist of 
scaling limits for the total spin and moderate deviation principles (MDPs) for the total spin. 

We will refer to the mean-field model studied in this paper as the BEG model; it is equivalent 
to the Blume-Emery-Griffiths model on the complete graph on n vertices. In contrast to the 
mean-field version of the Ising model known as the Curie-Weiss model, whose only phase 
transition is a continuous, second-order phase transition at the critical inverse temperature [17, 
§IV.4], the BEG model exhibits both a curve of continuous, second-order points; a curve of 
discontinuous, first-order points; and a tricritical point, which separates the two curves [22, 29]. 
It is one of the few models, and certainly one of the simplest, that exhibit this intricate phase- 
transition structure. 

Applications of the Blume-Emery-Griffiths model to a diverse range of physical systems are 
discussed in [22, §1] and in [29, §3.3], where the model is called the Blume-Emery-Griffiths- 
Rys model. As the latter reference points out, the model studied in the present paper is actually 
a mean-field version of a precursor of the Blume-Emery-Griffiths-Rys model due to Blume [3] 
and Capel [8, 9, 10]. With apologies to these authors, we follow the nomenclature of our earlier 
paper [22] by referring to this mean-field version as the BEG model. 
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The BEG model is defined by a probability distribution Pn,i3,K, where n equals the number 

of spins, P is the inverse temperature, and K is the interaction strength. We investigate the 
complex structure of the phase transitions in the model by deriving 36 different limit results for 
the total spin Sn as K) converges along appropriate sequences Kn) to points belonging 
to three separate classes: (1) the tricritical point, (2) the curve of second-order points, and (3) 
the single-phase region lying under that curve. In case 1, we obtain 13 scaling limits and 13 
MDPs; in case 2, 4 scaling limits and 4 MDPs; and in case 3, 1 scaling limit and 1 MDP. As 
we will see, the numbers 13, 4, and 1 represent natural and exhaustive enumerations of three 
classes of polynomials that arise in the related settings of the scaling limits and the MDPs. 

The existence of 18 = 13 + 4 + 1 scaling limits and 18 MDPs reflects the intricate structure 
of the phase transitions in the BEG model. It is hoped that our insights can also be applied 
to other statistical mechanical models that exhibit other types of phase transitions and critical 
phenomena and thus, presumably, other possibilities for scaling limits of macroscopic random 
variables like the total spin in the BEG model [19]. 

Before saying more about the limit theorems in the BEG model and their critical behavior, 
we summarize a number of facts concerning the phase-transition structure of the model [22]. 
For > and K > we denote by S^^k the set of equilibrium macrostates of the model 
corresponding to the macroscopic variable of the spin per site. In [22] it is proved that there 
exists a critical inverse temperature f3c = log 4 and that for /? > there exists a critical value 
Kc{(3) > having the following properties. 

1. For (3 > and < K < Kc{l3), S/s^k consists of the unique pure phase 0. 

2. For (3 > and K > Kc{(3), Ej3,K consists of two distinct, nonzero phases. 

3. For < /3 < /3c, as increases through Kc{P), Ep^x undergoes a continuous bifurcation, 
which corresponds to a second-order phase transition. 

4. For (3 > Pc, as K increases through Kc{P), S^^k undergoes a discontinuous bifurcation, 
which corresponds to a first-order phase transition. 

5. The point (/3c, Kc{/3c)) = (log 4, 3/[21og4]) in the positive quadrant of the (3-K plane 
separates the second-order phase transition noted in item 2 from the first-order phase 
transition noted in item 4. The point (/3c, Kc{Pc)) is called the tricritical point. 

The limit theorems to be considered in the present paper focus on the values of (3 and K 
in items 1, 3, and 5. For each such (/3, K), S^^k consists of the unique pure phase 0. Figure 1 
shows the corresponding portion of the phase diagram, which exhibits three sets A, B, and C. 
C is the singleton set containing the tricritical point {Pc, Kc{Pc)), B is the curve of second-order 
points defined by 

B = {(/3, K)eR^:Q<P<P,,K^ K,{P)}, (1.1) 
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and A is the single-phase region lying under B U C and defined by 

A = {{(3,K) eR^ -.0 < (3 < (3^,0 < K < K^{(3)}. 



(1.2) 




Figure 1: The sets A, B, and C 



In the remainder of this introduction we focus on the scaling limits for the total spin Sn when 
{P,K) converges to the tricritical point (/3c, -ft'c(/?c)) along appropriate sequences {l3n-,Kn). 
These scaling limits describe the limiting distribution of Sn/n^^^ with respect to Pn,^ri,Kn for 
appropriate choices of 7 G (0, 1/2). The simplest sequences for which the full range of scaling 
limits appear are defined in terms of parameters a > 0, > 0, 6 7^ 0, and A; 7^ by 



Pn = log(e^^ - 6/n") and = i^(/3„) - k/n' 



(1.3) 



where K{(3) = (e^ + 2)/(4/3) for (3 > 0. K{(3) coincides with K^{(3) for < /3 < /^^ and 
satisfies K{[3) > K^{[3) for P > 13^ [22, Thms. 3.6, 3.8]. A detailed overview of all the limit 
theorems in the paper, including those discussed here, is given in the next section. 

In each of the scaling limits the form of the limiting density reflects the influence of one or 
more of the sets A, B, and C that lie in a neighborhood of the tricritical point. The influence of 
those sets, which depends only on a and 9 and not on 6 or A; in (1.3), is shown in Figure 2. In 
that figure the positive quadrant of the a-9 plane is partitioned into the following sets. 



1. Three open sets labeled A, B, and C. 
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2. Three line segments labeled A + B, A + C, and B + C that separate the three open sets 
in item 1 . 

3. The point equal to (1/3, 2/3) and labeled A + B + C at which the three line segments in 
item 2 meet. 
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Figure 2: Influence of C, B, and A when (/?„, K„) (/3c, KdPc)) 

Figure 2 is a limit-theorem phase diagram that summarizes the critical behavior of the scal- 
ing limits in a neighborhood of the tricritical point. This critical behavior consists of the fol- 
lowing phenomena, which can be verified by examining the statement of the scaling limits in 
Theorem 7.1. 



1. When (a, 9) lies in one of the open sets labeled A, B, or C, then the limiting density in 
the corresponding scaling limit shows the influence only of that single set. Hence these 
three open sets correspond to the pure phases of the scaling limits. 

2. When (a, 9) lies in one of the line segments labeled A + B, A + C, or B + C, then 
the limiting density shows the influence of both sets, A and B, A and C, or B and C, 
respectively. Hence these three line segments correspond to the coexistence of the pure- 
phase scaling limits noted in item 1 . 

3. When {a, 9) equals the point labeled A + B + C, then the limiting density shows the influ- 
ence of all three sets A, B, and C . This point is the analogue of the tricritical point in the 
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standard phase diagram, a portion of which is shown in Figure 1. Indeed, any neighbor- 
hood of the tricritical point in the P-K plane contains values of (3 and K corresponding to 
all the different phase-transition behaviors of the model. Similarly, any neighborhood of 
the analogue of the tricritical point in the limit-theorem phase diagram contains values of 
a and 9 corresponding to all the different forms of the scaling limits, which number 13. 

4. As (a, 6) crosses any of the line segments labeled A + B,A + C,oxB + C, the values of 
7 in the scaling limits change continuously, which corresponds to a second-order phase 
transition; by contrast, the forms of the limiting densities change discontinuously, which 
corresponds to a first-order phase transition. 

As noted in items 1, 2, and 3, the influence of the sets upon the forms of the limiting den- 
sities reveals a fascinating geometric feature of the BEG model. This feature is completely 
unexpected because the model has no geometric structure. In fact, each spin interacts equally 
with all the other spins via a mean-field Hamiltonian, and so the model is independent of di- 
mension. The discussion of the scaling limits given here, including the notion of the influence 
of a set on the form of the limiting density, will be greatly amplified in the next section. 

The scaling limits of Sn/ n}''^ corresponding to the choices of a and 9 in Figure 2 are derived 
in Theorem 7.1, where we determine the values of a, 9, and 7 leading to the various forms of 
the limit. In Figure 2 the value or range of values of 7 are also shown for (a, 6) lying in the 
sets labeled A, B, and C. The set labeled A is divided into two subsets by the line 9 — 1/2; the 
ranges of values of 7 are different in the two subsets. 

The three seeds from which the present paper grew are [22], [20], and [14]. In the first paper 
the phase-transition structure of the BEG model is analyzed. In the second paper scaling limits 
are proved for a class of models that includes the Curie-Weiss model as a special case. In the 
third paper 4 different MDPs are obtained for the Curie-Weiss model when the inverse temper- 
ature converges to the critical inverse temperature in the model along appropriate sequences 
The results derived in the present paper greatly extend both the scaling limits in [20] and the 
MDPs in [14]. This is the case because the BEG model has a much more intricate structure of 
phase transitions than the Curie-Weiss model and so exhibits a much richer class both of scaling 
limits and of MDPs. As we will outline near the end of the next section, both the scaling limits 
and the MDPs are proved by a unified method. 

This unified method is based, in part, on properties of a function G^^k defined in (3.4). This 
function plays a central role in every aspect of the analysis of the BEG model considered in the 
present paper as well as in its prequel [22]. In summary these are the following. 

• The set Sp^K of equilibrium macrostates for the BEG model is defined as the set of zeroes 
of the rate function in the LDP for the P„ /3 /^-distributions of Sn/n given in Theorem 
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3.1. In turn, this set coincides with the set of global minimum points of Gp^K [see (3.5)]. 
This characterization of Epj^ allowed us to carry out the detailed analysis of the phase- 
transition structure of the model in [22]. 

• The canonical free energy K) equals the global minimum value of G^^k [see item 2 
after (3.4)]. 

• The distribution of Sn/v}~'^ can be expressed directly in terms of G^^k [Lem. 4.1]. 

• Gjj^K is the rate function in a second LDP involving Sn/n given in part (b) of Lemma 
4.4. The estimates derived from this LDP and given in parts (c) and (d) of the lemma are 
the key estimates needed to control error terms in the proofs of the scaling limits and the 
MDPs for Sn/n}"^- Lemma 4.4 is the main technical innovation in the paper. 

• When a certain quantity w defined in terms of a, 6, and 7 equals 0, the 13 different 
forms of the Taylor expansion of nG p^^K„{x / rf) for appropriate sequences Kn) and 
7 e (0, 1/2) yield the 13 different forms of the scaling limits of Sn/n}''^ [Thm. 7.1]. 

• When w < 0, the 13 different forms of the Taylor expansion of 'n}'^'^Gf3^^K„{x/n'^) for 
appropriate sequences Kn) aad^ e (0, 1/2) yield the 13 different forms of the MDPs 
of Sjn^-^ [Thm. 8.3]. 

This discussion shows that all the magic is in the function Gp^x- The fact that the wide 
variety of phenomena derived in the present paper and in [22] can be obtained via properties of 
a single function is an appealing feature of the BEG model. Besides the Curie- Weiss model and 
generalizations studied in [14, 20, 21, 30] and numerous other papers, this feature is shared with 
a number of other mean-field models, including a mean-field version of the nearest neighbor 
Potts model known as the Curie-Weiss-Potts model [25], the mean-field XY Heisenberg model 
[1], and the Hopfield model of spin glasses and neural networks [31]. These mean-field models 
have in common the fact that the interaction terms in their Hamiltonians can be written as a 
quadratic function. Scaling limits and MDPs for these models have either been proved, or 
in principle could be proved, by techniques similar to those used in the present paper. Some 
of these techniques are generalized in [11], in which the quadratic term in the Hamiltonian is 
replaced by the moment generating function of suitable random variables. Other generalizations 
are given in [5, 6, 23, 24]. The analysis of the equilibrium macrostates and the associated 
phase transitions in the BEG model, which underlies the present paper, is carried out in [22] 
using large deviation techniques. While this is an elegant method that provides exact, analytical 
results, it has the restriction that it works most efficiently in models with long-range interactions, 
as explained in [2]. 
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The Hopfield model of spin glasses and neural networks has received a great deal of at- 
tention, and limit theorems for this model have been actively studied. The Hamiltonian in the 
Hopfield model can be written as a quadratic function of the overlap parameter, a feature that 
it shares with the Curie-Weiss model and the BEG model, in which the Hamiltonian can be 
written as a quadratic function of the spin per site. For the Hopfield model both central limit 
theorems and non-classical scaling limits for the overlap parameter are studied in [7, 26, 27, 28], 
and MDPs are studied in [15]. These limit theorems include the cases when the inverse temper- 
ature is constant and when the inverse temperature parameter converges to the critical inverse 
temperature at an appropriate rate [15, 28]. 

We next preview the contents of the present paper. In section 2 a detailed overview is 
given of the scaling limits and the MDPs that will be derived. In section 3 we summarize the 
results in [22] on the structure of the set of equilibrium macrostates of the BEG model and the 
associated phase transitions. In section 4 we introduce the function Gp^x, properties of which 
are integral to the proofs of the scaling limits and MDPs. These properties include a formula 
for the distribution of the total spin in terms of Gp^K [Lem. 4.1], several forms of the Taylor 
expansions of Gp^K that will be used to derive the limit theorems [Thm. 4.3], and two estimates 
in Lemma 4.4 for controlling error terms in the proofs of the scaling limits and the MDPs. 

In sections 5-8 we apply the results in the previous sections to derive the scaling limits 
and the MDPs. Sections 5 and 6 are devoted to scaling limits for Sn/n}~^ when appropriate 
sequences (Pn^Kn) converge to points {(3,K) e A and to points {/3,Kc{(3)) e B, where A and 
B are the sets defined in (1.2) and (1.1). When Kn) {P, K) e Awe obtain only 1 scaling 
limit, which is independent of the sequence Kn) [Thm. 5.1]. The situation for (/?, Kc{(3)) G 
B is much more interesting; for appropriate choices of Kn) (P, K) G B, 4 different 
forms of the scaling limits arise [Thm. 7.1]. The scaling limits proved in these two sections 
are warm-ups for the even more complicated scaling limits proved in section 7. In that section, 
for appropriate choices of Kn) converging to the tricritical point (/3c, Kc{Pc)) we obtain 
13 different forms of the scaling limits [Thm. 7.1]. Finally, in section 8 we obtain 1 MDP for 
Sn/n^-^ when (/3„, Kn) (P, K) e A [Thm. 8.2], 4 MDPs when (/3„, Kn) ^ (/3, K,{(3)) e 
B [Thm. 8.1], and 13 MDPs when (/3„, Kn) {(3c, Kc{(3c)) [Thm. 8.3]. The MDPs are proved 
by showing the equivalent Laplace principles, which is carried out by a method closely related 
to that used to prove the scaling limits in the earlier sections. Being able to prove both classes 
of limit theorems via a unified method is one of the attractive features of this paper. 

Acknowledgements. We would like to thank an anonymous referee of [22] who suggested 
studying scaling limits for Sn/n}~'^ in the BEG model using sequences (/3„, Kn) converging to 
various points (/3, K). We would also like to thank Jonathan Machta for useful discussions on 
the material of the present paper. The research of Richard S. Ellis is supported in part by a grant 
from the National Science Foundation (NSF-DMS -0604071). The research of Peter Otto was 
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2 Overview of the Limit Theorems 

This paper is devoted to scaling limits and MDPs for the total spin in the BEG model. In 
order to highlight the novelty of these results, we introduce some notation. The BEG model 
is a lattice-spin model defined on the complete graph on n vertices 1, 2, . . . , n. The spin at 
site j e {1, 2, n} is denoted by loj, a quantity taking values in A = {—1, 0, 1}. The joint 
distribution of the spins LOj is defined by a probability measure Pn,0,K on the configuration 
space A" [see (3.1)]. The sequence Pn,f3,K forn e N defines the canonical ensemble for the 
BEG model. 

Through the particular form of the interactions among the spins, the measures Pn,i3,K incor- 
porate an alignment effect that underlies the phase-transition structure of the model. As /3 — > 0, 
Pn,i3,K converges weakly to the product measure on A" with marginals equal to the uniform 
measure on A. Similarly, as 7^ ^ 0, Pn.g.K converges weakly to another product measure on 
A". By contrast, as A' ^ oo, Pns,K concentrates on the configurations uj^ and uj^ in which the 
spins are all 1 or —1; by symmetry, as AT — cxd, Pn,p,K converges weakly to the sum of point 
masses | {6^+ + 5^-). The phase-transition structure of the model reflects the persistence of this 
alignment effect in the limit n — > oo. 

We define Sn = S"=i <^j^ which represents the total spin. In this paper we will consider 
numerous weak limits of the distributions of Sn/n^~"', where 7 G [0, 1). The distributions are 
with respect to Pn,f3,K for fixed /3 > and AT > and, more generally, with respect to Pn,f3„,Kn'> 
where (/3„, A'„) are appropriate sequences converging to specific values of {/3, K). The use of 
Pn,i3n,K„ to Study weak limits in place of Pn,f3,K is the basic innovation of this paper, which will 
reveal the intricate phase-transition structure of the model. If z/ is a probability measure on M, 
then the notation Pn,(3„,K„{Sn/n'^^^ ^ dx} =^ v means that the distributions of Sn/n^^^' with 
respect to Pn,0n,K„ converge weakly to z/ as n — > cxo. If / is a nonnegative integrable function 
on M, then the notation Pn,i3.a,Kn{Sn/n^~'^ G dx} =^ f dx means that the distributions of 
Sn/n^~'^ converge weakly to the probability measure on M having a density proportional to / 
with respect to Lebesgue measure. 

The first hint of the intricacy of the phase-transition structure of the BEG model can be seen 
by examining the law of large numbers and its breakdown, which we consider with respect to 
Pn,/3,K for fixed /3 > and K > 0. The intuition is that for sufficiently small X > the 
interactions among the spins are sufficiently weak so that the analogue of the classical law of 
large numbers holds. However, for sufficiently large K > the interactions among the spins 
are sufficiently strong to cause the classical law of large numbers to break down. This intuition 
is in fact correct. In [22] it is proved that there exist Kc{(3) > 0, defined for /3 > Q, and z{l3, K), 
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defined for /3 > and K > Kc{0), in terms of which the following limits hold. The form of 
the limits for K = Kc{(3) is given in (2.3) and (2.4). 

• For any /3 > and < X < K^{I3) 

Pn,i3,K{Sn/n G dx} =^ ^o- (2.1) 

• For any /3 > and K > K^P) we have z{P, K) > and 

Pn,/3,K{Sn/n G dx} l{Sz{f3,K) + S-z{i3,K)) ■ (2.2) 

The proofs of these two limits are indicated at the end of section 3, where they are derived from 
the LDP given in part (a) of Theorem 3.1. 

As we explain in section 3, for each P > and > the sets of mass points of the limiting 
measures represent the sets of equilibrium macrostates of the BEG model, which we denote by 

Sp^K- Thus, for /5 > and < < K^{f3), Sp^x = {0} while for /3 > and > K^{p), 
^(3,K = {^z{f3, K)}. The quantity z{l3,K) is a positive, increasing, continuous function for 
K > K^{I3). The limit of z{j3, K) as K ^ K^{I3)+ depends on whether p < f3c or p > p^, 
where Pc — log 4 represents the critical inverse temperature of the model. For P > P^we have 
z{P, K,{P)) > 0, and 



K 



^^ _ / ifO < /3 < Pc 

" 1 zip, K^m ifp> Pc. 



Consistent with this limit behavior is the fact that Sp^Kd/s) equals {0} for < /3 < /3c and equals 
{0, ±z{p, Kc{P))} for p > Pc. The limit behavior of z{p, K) exhibited in the last display 
shows that the sets £p^K undergo a continuous bifurcation dX K — Kc{P) for < P < Pc 
and a discontinuous bifurcation at = K^P) for P > Pc. From the viewpoint of statistical 
mechanics, the dual bifurcation behavior of the model corresponds to a continuous, second- 
order phase transition at {P,Kc{P)) for < P < Pc and a discontinuous, first-order phase 
transition at {p, KdP)) for p > Pc. The point {Pc, K^pc)) = (log 4, 3/[21og4]) separates the 
second-order phase transition from the first-order phase transition and is called the tricritical 
point. 

The different behavior of the two phase transitions is reflected in the form of the limits of 
Sn/n when K = Kc{P). For < /5 < /5c, we have the law of large numbers 

PnAKc{P){Sn/n e dx} =^ 5o, (2.3) 

while for P > Pc the limit is expressed in terms of a measure supported at the three points in 

Pn,0,Kc{(3){Sn/n G dx} =^ \q5q + \l{5^^p^K,(0)) + (^-2(/3,Kc(/9))) ■ (2-4) 
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In the last limit Aq and Ai are positive numbers satisfying Aq + 2Ai = 1 and given explicitly in 

(4.4). As we point out at the end of section 3, (2.3) follows immediately from the LDP given in 
part (a) of Theorem 3.1. However, the proof of (2.4) is more subtle and is postponed until after 
Theorem 4.2. 

Further evidence of the intricacy of the phase-transition structure of the model can be seen 
if one jumps from the context of the law of large numbers and its breakdown to the context of 

scaling limits for Sn that are related to the central limit theorem and its breakdown. We consider 
three cases, in all of which S^^k = {0}. Case 1 is defined by /? > and < K < Kc{P). For 
these values of (3 and K the interactions among the spins are sufficiently weak, and the analogue 
of the classical central limit theorem holds. As we prove in Theorem 5.1 when < /3 < (3c, 

Pn,i3,K{Sn/n^^'^ G dx} =^ exp(-C2X^) dx, (2.5) 

where C2 = C2{P,K) is defined in (5.1). The same limit holds when P > Pc and < K < 

Kci(3). 

Case 2 is defined hy < P < Pc and K = Kc{P). In this case the central limit scaling n^/^ 
in (2.5) must be replaced by n^^*, which reflects the onset of long-range order represented by 
the second-order phase transition at {P, Kc{P)). We have the nonclassical limit 

Pn,i3,Kc{l3){Sn/n^^'^ G dx} =^ exp(-C4X^) dx, (2.6) 

where C4 = c^iP, K) > is defined in (6.5). The limit in the last display is a special case of 
one of the limits proved in Theorem 6. 1 [see the note after the statement of the theorem]. 

Case 3 focuses on the tricritical point (/9c, Kc{Pc))- Not only is there an onset of long-range 
order represented by the second-order phase transition at this point, but also this point separates 
the second-order phase transition for /3 < /3c and the first-order phase transition for /? > Pc. This 
more intricate phase-transition behavior in a neighborhood of the tricritical point is reflected in 
the replacement of the scaling n^/^ for < /9 < /9c by rv'/^. In this case 

Pn,i3MI3o){Sn/n^'^ £ dx] =^ exp(-C6X^) dx, (2.7) 

where ce = 9/40. The limit in the last display is a special case of one of the limits proved in 
Theorem 7.1 [see the note after the statement of the theorem]. 

For all other values of P > and K > — those satisfying < P < P^., K > Kc{P) and 
P > Pc, K > Kc{P) — the limit theorems have different forms because the set Ep^K of equi- 
librium macrostates consists of more than one point. In both of these cases, for any equilibrium 
macrostate z, {Sn — nS)/n^/^ satisfies a central-limit-type limit when Sn/n is conditioned to lie 
in a sufficiently small neighborhood of z. The explicit form of the limit is given in part (b) of 
Theorem 6.6 in [22]. 
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We are now ready to outline the main contribution of this paper, which is to exhibit the 

intricate probabilistic behavior of the BEG model in neighborhoods of the tricritical point 
{l3c, Kc{l3c)), second-order points (/?, Kc{l3)) for < /5 < /5c, and points {/3, K) for < f3 < 13^ 
and Q < K < Kc{(3). We do this by studying scaling limits and MDPs for Sn/in}~^ with re- 
spect to Pn,/3„,K„ for appropriate sequences K^) that converge to points belonging to these 
three classes and for appropriate choices of 7 e (0, |]. In order to facilitate the discussion, we 
denote by C the singleton set containing the tricritical point (/3c, Kc{(5^), by B the curve of 
second-order points defined by 

B = {{(5,K) e . < /3 < /3e,ii: = i^c(/3)}, 

and by A the single-phase region lying under B\JC and defined by 

A = {{(5,K) e < /3 < /3e,0 < X < K,{(5)). 

The sets A, B, and C are shown in Figure 1 in the introduction. In the rest of this section 
we focus mainly on the scaling limits and MDPs for Sn/'n}~'^ when (/3„, K^) is an appropriate 
sequence that converges to {/3c, Kc{Pc))- Scaling limits and MDPs when (/5„, Kn) converges to 
{(3, Kc{(3)) E B and to (/?, K) E A dXQ treated, respectively, in Theorems 6.1 and 8.1 and in 
Theorems 5.1 and 8.2. 

Corresponding to each {(3, K) e AliBliC there exists a unique equilibrium macrostate at 0. 
We do not consider scaling limits and MDPs in the neighborhoods of other points corresponding 
to which there exist nonunique equilibrium macrostates. In all or most cases of nonunique 
equilibrium macrostates, we expect that the scaling limits and MDPs are conditioned limits as 
in [22, Thm. 6.6(b)] and [14, Thm. 1.1]; however, we have not worked out the details. 

Through the limits (2.5), (2.6), and (2.7), each of the sets A, B, and C is associated, respec- 
tively, with the term x^, x^, and x^. Specifically, for fixed (/3, K) 

{exp{-C2X^) dx with-f= 1/2 if {f3,K) e A 
exp(-C4a;4) dx with 7 = 1/4 if (/3, K) e B (2.8) 
exp(-C6x6) dx with 7 = 1/6 if (/3, K) e C, 

where C2 and C4 are positive and depend on P and K, and cg = 9/40. Theorem 7.1 shows 
that for appropriate sequences {(3n-,Kn) converging to (/3c, -f^c(/3c)), for appropriate choices of 
7 e (0, 1/2), and for appropriate coefficients 62, 64, and ce 

Pn,0n,Kn{Sn/'n^~'^ e dx} =^ exp{-C2x'^ - c^x'^ - cqx^) dx. (2.9) 

As we show in Table 2.1, G{x) = C2X^ + + c^x^ takes all of the 13 possible forms of an 
even polynomial of degree 6, 4, or 2 satisfying G(0) — and G{x) ^ 00 as \x\ 00. Each of 
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the 13 cases shows the influence of one or more of the sets C, B, and A through the presence 
of the term x^, x^, or x"^ associated with that set by the limit (2.8). The coefficient cg = 9/40 is 
the same as in (2.7), q = 3/16, and b and k are any nonzero real numbers subject only to the 
requirement that exp(— G) is integrable. Because in every case 7 G (0, 1/2), the scaling of Sn 
by n^"'"' is non-classical. 



case 


influence 


Prv,f3„,K„{S.n/n^ ^ G dx} ^> exp[-G{x)] dx 


1 


C 


G{x) = CqX^, Cq > 


2 


B 


G{x) = hcix'^, 6 > 0, C4 > 


3 


A 


G{x) = k(3cx'^, k>0 


4-5 


B + C 


G{x) = bc4X^ + C6X^ by^O 


6-7 


A + C 


G{x) = kAx'^ + cgx^, k^O 


8-9 


A + B 


G{x) = k/3cx^ + bc4x\ k^O,b>0 


10-13 


A + B + C 


G{x) = kPcX^ + bc^x"^ + cgx^ k^O,b^O 



Table 2.1: 13 cases of the scaling limits in (2.9) for Kn) in (2.10) and 7 G (0, 1/2) 



The forms of the scaling limits in (2.9) depend crucially on the appropriate choices of the 
sequences (/5„, Kn) converging to (/5c, Kc{Pc))- The simplest sequences for which all 13 cases 
of the limit (2.9) arise are defined in terms of parameters a > 0, ^ > 0, 6 7^ 0, and k 

Pn = logie"^ - b/n'') and = K{Pn) - k/n\ (2.10) 

where K{p) = (e^ + 2)/(4/3) for (3 > 0. K{(3) coincides with K^((3) for < < /^^ and 
satisfies K{P) > K^iP) for p > pc [22, Thms. 3.6, 3.8]. Since pn pc and since K{-) is 
continuous, we have K{Pn) K^, thus the convergence (/9„, Kn) {Pc, Kc{Pc)) is valid. In 
section 7 we will explain how this particular sequence {Pn,Kn) was chosen. 

Depending on the signs of b and k, the sequence (/3„, Kn) in (2.10) converges to {P^ Kc{Pc)) 
from regions exhibiting markedly different physical behavior. For example, if 6 > and k > 
0, then Pn < Pc and Kn < K{Pn), and so (Pn, Kn) converges to iP.K) from the region 
A, corresponding to each point of which there exists a unique equilibrium macrostate [Thm. 
3.2(a)]. On the other hand, if /c < 0, then Kn > K{Pn), and so {Pn, Kn) converges to {P, K) 
from a region of points corresponding to each of which there exist two equilibrium macrostates. 
If, in addition, 6 > 0, then this region lies above the curve B of second-order points [Thm. 
3.2(b)], while if 6 < 0, then this region lies above the curve of first-order points described in 
Theorem 3.3. Despite the markedly different physical behavior associated with these various 
regions, all the scaling limits in this paper are proved by a unified method, regardless of the 
direction of approach of {Pn, Kn) to {P,K). The situation with respect to the MDPs is the same. 
These remarks concerning the proofs of the scaling limits and the MDPs will be amplified in 
section 4 after we introduce the tools that will be used in the proofs. 
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The occurrence of a particular one of the scaling limits enumerated in Table 2. 1 depends on 

7 and on the values of a. and 9 and thus on the speed at which Kn) {0c, Kc{Pc)) and on 
the direction of approach. Only case 1 expresses the influence of C alone, giving the same limit 
for (Pn, Kji) in (2.10) as the limit in (2.7), which holds for the constant sequence Kn) — 
{Pc: ^ciPc))- Case 1 occurs if the convergence (/3„,ir„) — > {Pc: ^^dPc)) is sufficiently fast; 
namely, a > 1/3 and 6 > 2/3. Case 2, which expresses the influence of B alone, occurs if the 
convergence of Kn) {Pc, -^c(/5c))is sufficiently slow but 9 is relatively large compared 
to a.. Case 3, which expresses the influence of A alone, occurs if the convergence is sufficiently 
slow but, in contrast with case 2, a is relatively large compared to 9. Finally, cases 4-13, 
which express the influence of more than one set A, B, and C, occur if the convergence of 
{Pn: Kn) — > (/3c, Kc{Pc)) occurs at an appropriate critical rate. For example, cases 10-13 
express the influence of all three sets A, B, and C and so correspond to the most complicated 
form of the limiting density. This case occurs when a = 1/3,9 = 2/3, and 7 = 1/6. 

The scaling limits for Sn/n^~^ listed in Table 2.1 are analyzed in Theorem 7.1, where we 
determine the values of a, 9, and 7 leading to the 13 different cases. The dependence of K^) 
in (2.10) upon a and 9 is complicated; because Pn is a function of a, Kn is both a function of 
9 and, through Pn, a function of a. However, as we will see, for the appropriate choice of 
7 G (0, 1/2), in the expression for the scaling limit of Sn/n^^^' the a and the 9 decouple in such 
a way that the limits given in Theorem 7. 1 can be read off in a systematic way. 

In Figure 2 in the introduction we indicate the subsets of the positive quadrant of the a-9 
plane leading to all the cases in Table 2.1. The subsets labeled C, B, and A correspond to 
cases 1, 2, and 3, respectively, and the subsets labeled B + C, A + C, A + B, and A + B + C 
correspond to cases 4-5, 6-7, 8-9, and 10-13, respectively. The relationship between the a-9 
plane exhibited in Figure 2 and the P-K plane, inside which lies the tricritical point, is that each 
point in the a-9 plane corresponds, through the formulas for /3„ and Kn given in (2.10), to a 
curve in the P-K plane. 

In Figure 3 we exhibit three different curves in the P-K plane, labeled (a), (b), and (abc). 
These curves correspond to three different choices of a and 9, three different choices of (Pn, K^) 
in (2.10), and three different limits in Table 2.1. The curve labeled (a) corresponds to a = 1 
and 9 — 1/3, which in turn corresponds to case 3 of the scaling limit; this case shows the 
influence only of region A. The curve labeled (b) corresponds to a = 1/4 and 9 — 1, which 
in turn corresponds to case 2 of the scaling limit; this case shows the influence only of region 
B. Finally, the curve labeled (abc) corresponds to a = 1/3, 6 > 0, = 2/3, and A; > 0; the 
associated scaling limit in case 10 shows the influence of all three sets A, B, and C. 

It is worth noting a contrast between the scaling limits in (2.8) and those in Table 2.1. In 
(2.8) the three scaling limits for Sn/'n}~'^ hold with respect to Pn,i3,K for fixed K) e A, 
{P, K) E B, and (/?, K) G C. In each of these three cases the value of 7 is fixed to be, 
respectively, 1/2, 1/4, and 1/6. By contrast, we will see in Theorem 7.1 that in 4 of the 13 
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Figure 3: Three choices of K„) that show the influence of A, of B, and of A, B, and C in (2.9) 



cases of the scaling limits for Sn/n^~'^ stated in Table 2.1, the limit theorems hold for a range 
of values of 7. These are cases 2, 3, 8, and 9. In the other cases, each of which includes the 
influence of the tricritical point (/9c, Kc{Pc)), 7 equals the fixed value 1/6. 

We now make a transition from the scaling limits to the MDPs. As we have seen, the scaling 
limits state that for appropriate choices of (/3„, Kn) and of 7 = 70 

PnAuKASn/n'-'"' e dx} =^ exp[-G'(x)] dx, (2.11) 

where G takes one of the 13 forms in Table 2.1. For any 7 G (0, 70), one can show that if D is 
any Borel set whose closure does not contain 0, then 

lim Pn,f3„,KASn/n'-'' E D} = 0. 

n— >oo 

A natural question is to determine the rate at which these and related probabilities converge to 
when Kn) is defined in (2.10). In Theorem 8.3 we define a quantity w in terms of a, 
9, and 7 having the property that when w < 0, Sn/n^^'^ satisfies an MDP with exponential 
speed n~'^ and rate function G{x) — G, where G is the same function appearing in (2.11) and 
G = infj^gK G{y). This MDP implies that for suitable sets D 

Pn,p^,KASnln^-^ G ^ like expl-n-"- inf - G)]. 

In order to emphasize the similarity with the scaling limits, we summarize this class of MDPs 
by the formal notation 

Pn,p^,KASnn'~^' G dx} X exp[-n--G(x)] , (2.12) 
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in which the constant G is not shown. 

The situation with the MDPs is completely analogous to the situation for the scaling limits. 
Specifically, as we exhibit in Table 2.2, there are 13 cases of the MDP (2.12), each of which 
shows the influence of one or more of the sets C, B, and A depending on the speed at which the 
sequence Kn) defined in (2.10) converges to (/3c, Kc{Pc)) and on its direction of approach. 
The coefficient cq — 9/40 is the same as in (2.7), C4 = 3/16, and b and k are the nonzero real 
numbers appearing in (2.10) and subject only to the requirement that G{x) ^ 00 as 00. 
The MDPs for Sn/n^^"' listed in Table 2.2 are analyzed in Theorem 8.3, where we determine 
the values of a, 9, and 7 that lead to each of the cases. 



case 


influence 


Pr^,/3„,K„{5„/nl-^ ^ ^ exp[-n— G(a;)] dx 


1 


G 


G{x) = CqX^, Cq > 


2 


B 


G{x) = hc^x^, 6 > 0, C4 > 


3 


A 


G{x) = kl3,x^ k>0 


4-5 


B + C 


G{.r) = hcix^ + (:(iX^\ h ^ 


6-7 


A + G 


G(x) = kPcX^ + cqx®, /c 7^ 


8-9 


A + B 


G{x) = k(3cx^ + 6c4x\ A; ^ 0, 6 > 


10-13 


A + B + G 


G{x) = kl3cX2 + hcis^x* + Cqx^, k 0, b 



Table 2.2: 13 cases of the MDPs in (2.12) for (/3„, Kn) in (2.10) and 7 G (0, 1/2) 



The MDPs for Sn/n^^'^ have an unexpected consequence concerning a new class of distri- 
bution limits for Sn/n^~'^ that give deeper insight into the fine structure of the phase transitions 
in a neighborhood of the tricritical point. In an effort to understand the physical significance of 
these new limits, analogs of them are now being investigated for a class of non-mean-field mod- 
els, including the Blume-Emery-Griffiths model [19]. In order to appreciate these new results, 
we first consider a consequence of the large deviation principle stated in part (a) of Theorem 
3.1. Since Sp^K = {0} for {(3,K) e A U B U G, it follows that for any positive sequence 
{/3n,Kn)^{P,K)eAuBUG 

Pn,f3„,Kn{Sn/n G dx} =^ 60- 

The MDPs for Sn/n^~"' listed in Table 2.2 lead to refinements of this limit for (Pc, Kc{(3c)) E 
G in those cases in which the set of global minimum points of G contains nonzero points. These 
are precisely the cases in which the coefficients of G are not all positive: cases 5 (6 < 0), 7 
(k < 0), 9 (A; < 0), 11 (A; < 0, 6 > 0), 12 (A; > 0, 6 < 0), and 13 (A; < 0, 6 < 0). In all these 
cases except for case 12, the set of global minimum points of G consists of two symmetric, 
nonzero points ±x(b, k). Hence, using the appropriate value of 7 and the appropriate sequence 
(/?„, Kn) given in Theorem 8.3, we deduce from the corresponding MDP the limit 

-Pn,/3„,K„{'S'„/n^"^ e dx] =^ \ {5^{b,k) + (^-x(6,fc))- (2.13) 
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In each of these cases {Pn,Kn) approaches {Pc, Kc{Pc)) from a region of points {l3,K) cor- 
responding to each of which there exist two equilibrium macrostates ±z{(3, k) [Thms. 3.2(b), 
3.3(c)]. As we have already seen, for each {P,K) in this region the limit (2.2) holds. The 
new limit (2.13) shows that as Kn) {/3c, Kc{(3c)) from this two-phase region, the model 
retains a trace of the two equilibrium macrostates ±^(/3, K), replacing them by the quantities 
±x{b, k). The physical significance of this limit as well as the limit (2.14) to be stated in the 
next paragraph is currently under investigation [19]. A similar phenomenon occurs in case 4 of 
Theorem 8.1, which proves MDPs for Sn/n^^'^ for appropriate sequences Kn) converging 
to {13, K) lying in the curve B of second-order points. 

The situation in case 12 in Table 2.2 (A; > 0, 6 < 0) is even more fascinating than in the other 
cases. For fixed 6 < 0, fixed n e N, and decreasing A; > 0, the set of global minimum points 
of G undergoes a discontinuous bifurcation, changing from a unique global minimum point at 
for k large, to three global minimum points at 0, ±,x(6, k) for a critical value oik = const ■ 6^, 
to two global minimum points at ±x{h, k) for k small. As k decreases, Kn) crosses the 
first-order critical curve from below; the changing forms of the sets of global minimum points 
of G replicate the changing forms of Ejj^K for fixed P > Pc and increasing K > [Thm. 
3.3]. In particular, when the set of global minimum points of G equals {0, ±x{b, k)}, the MDP 
corresponding to case 12 together with other information yields the limit 

Pn,(3„,KASn/n^'~'^ G dx} ^ XqSq + >^l{Sx{b,k) + S-x(b,k)): (2.14) 

where Aq and Ai are positive numbers satisfying Aq -I- 2Ai = 1. This limit is reminiscent of the 
limit (2.4), in which the equilibrium macrostates ±z{P, K) are replaced by their traces ±x(6, k). 

Although in general the values of a, 9, and 7 leading to each of the 13 cases of the MDPs 
in Table 2.2 differ from the values of these parameters leading to the corresponding case of the 
scaling limit in Table 2.1, the tables have a number of obvious similarities. This resemblance 
between the two tables reaches deeper. In fact, both sets of results are proved by a unified 
method. In order to explain this, let / be any bounded, continuous function mapping M into M 
and let (/3„, Kn) be any positive sequence. The starting point of the proofs of both the scaling 
limits and the MDPs [see Lem. 4.1] is that whenever 7 e (0, 1/2), we have 

E{f{Sn/n'-'' + En)} = ^ • / f{x) exp[-nG^„,K„(x/n^)] dx. (2.15) 

The function G/s^k in this display is defined in (3.4); its global minimum value equals the 
canonical free energy for the model. In addition, represents a sequence of random variables 
that converges to as n — > 00, and Z„ is a normalizating constant. 

The quantity w in the MDP (2.12) is defined hy w = min{27 -|- ^ — 1, 47 + a — 1, 67 — 1}. 
This quantity also plays a key role in the scaling limits for Sn/n}''^, which like the MDPs arise 
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from the choice of Kn) in (2.10). When w — 0, the scaling limits listed in Table 2.1 follow 
at least formally from (2.15) and the fact that for each x e IR 

lim uGp^^kSx/v?) = G{x). 

n— >oo 

The proof of this limit relies on an analysis of the Taylor expansion of Gf3^^K„ at 0, which has 
13 different forms depending on the choices of 7 and of the parameters a and 9 appearing in 
the definition (2.10) of Kn). Details are given in Theorem 7.1. 

We now assume that w < 0. Given t/j be any bounded, continuous function, we substitute 
/ = exp(n~"''0) into (2.15), obtaining 

^ 1 . /■ exp[n-- {ij{x) - n'+^Gp^,KAx/n^)}] dx. 

When w < 0, the last display, the fact that for each x e R 

lim n'+'"G0^,K„(x/n'^)^G{x), 

n—foo 

and the fact that — > in probability at a rate faster than exp{—n~'^) give the formal asymp- 
totic s 

E{exp[n-'"ij{Sn/n'-^ + Sn)]} 

?a [ exp[n-"' {^Pix) - {G{x) - G)}] dx 

~ exp [n""' sup^gK {ij{x) - {G{x) -G)}], 

where G = infj^gK G{y). In section 8 we show to convert this formal calculation into a limit 
known as the Laplace principle, which is equivalent to the MDPs for Sn/n^~'^ listed in Table 
2.2. As in the proof of the scaling limits, the proof of the Laplace limit relies on an analysis of 
the Taylor expansion of Gp^^K„ at 0. Despite the similarity in the proofs of the scaling limits 
and the Laplace principles, the proof of the latter is much more delicate, requiring additional 
estimates not needed in the proof of the former. 

We start our analysis of the BEG model in the next section. 



3 Phase-Transition Structure of the BEG Model 

After defining the BEG model, we summarize its phase-transition structure in Theorems 3.2 
and 3.3. In (3.4) we introduce the function Gp^x, in terms of which the scaling limits and the 
MDPs for Sn/n}''^ will be deduced later in the paper. 
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The BEG model is a lattice-spin model defined on the complete graph on n vertices 1 , 2, . . . , n. 
The spin at site j E {1, 2, . . . , n} is denoted by cuj, a quantity taking values in A = {—1, 0, 1}. 
The configuration space for the model is the set A" containing all sequences cu = (ui, uj2, ■ ■ ■ , ujn) 
with each uj G A. In terms of a positive parameter K representing the interaction strength, the 
Hamiltonian is defined by 



-I 




for each cu e A". For n E N, inverse temperature (3 > 0, and K > 0, the canonical ensemble 
for the BEG model is the sequence of probability measures that assign to each subset B of A" 
the probability 

PnAK{B) = ^j^^ ■ ^ ^M-PHn,K] dP^. (3.1) 

In this formula P„ is the product measure on A"^ with identical one-dimensional marginals 

p= |(5_l + 5o + (5i), 

and Zn{P, K) is the normalizing constant ex'^[—(3Hn,K\dPn- 

In [22] the analysis of the canonical ensemble Pn,/3,K was facilitated by expressing it in the 
form of a Curie-Weiss-type model. This is done by absorbing the noninteracting component of 
the Hamiltonian into the product measure P„, obtaining 



Pn,f3,K{duj) = ~ ■ exp 



nf3K 



PnAd^)- (3.2) 



In this formula Sn{uj) equals the total spin Yl^i^j^ Pn,0 is the product measure on A" with 
identical one-dimensional marginals 

pp{dwj) = ■ exp(-/3w|) p{dujj), (3.3) 

Z{I3) is the normalizing constant exp{—Pujj)p{dujj) = 1 + 2e~^, and Zn{l3, K) is the nor- 
malizing constant [Z{l3)f / Z^{I3, K). 

Although Pn,(3,K has the form of a Curie-Weiss model when rewritten as in (3.2), it is much 
more complicated because of the /9-dependent product measure Pn,/3 and the presence of the 
parameter K. These complications introduce new features not present in the Curie- Weiss model 
[17, §IV.4, §V.9]; these features include the existence of a second-order phase transition for all 
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sufficiently small (3 > and all sufficiently large K > and a first-order phase transition for 

all sufficiently large /? > and all sufficiently large 7^ > 0. The existence of a second-order 
phase transition and a first-order phase transition also implies the existence of a tricritical point, 
which separates the two phase transitions and is one of the main focuses of the present paper. 

The starting point of the analysis of the phase-transition structure of the BEG model is the 
large deviation principle (LDP) satisfied by the spin per site Sn/n with respect to Pn,p,K- In 
order to state the form of the rate function, we introduce the cumulant generating function cp of 
the measure pp defined in (3.3); for t e M this function is defined by 



We also introduce the Legendre-Fenchel transform of Cf^, which is defined for 2; e [—1, 1] by 

Jp{z) = sup{tz - Cfs{t)}; 



Jfsiz) is finite for z G [—1, 1]. Jp is the rate function in Cramer's theorem, which is the LDP for 
Sn/n with respect to the product measures Pn^p [17, Thm. II.4.1] and is one of the components 
of the proof of the LDP for Sn/n with respect to Pn,i3,K- This LDP and a related limit are stated 
in parts (a) and (b) of the next theorem. Part (a) is proved in Theorem 3.3 in [22], and part (b) 
in Theorem 2.4 in [18]. 

Theorem 3.1. For all (3 > and K > the following conclusions hold. 

(a) With respect to the canonical ensemble Pn,/3,K, Sn/n satisfies the LDP on [—1, 1] with 
exponential speed n and rate function 



where Zn{P, K) is the normalizing constant in (3.1). Then (p{(3, K) = miyQ^{Jfj{y) — j3Ky'^}. 

The LDP in part (a) of the theorem implies that those z G [—1,1] satisfying Ip^Kiz) > 
have an exponentially small probability of being observed in the canonical ensemble. Hence 
we define the set of equilibrium macrostates by 





(b) We define the canonical free energy 

ip((3, K)^- lim 



\ogZn{(5,K), 



£p,K^{ze[-l,l]:Ip,K{z)^Q}. 
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In [22] we used the notation S^^k to describe this set, using the notation Sp^x to describe a 
different but related set of equilibrium macrostates. In the present paper we write Sp^x instead 
of £f3^K in order to simplify the notation. 
For 2; e M we define 

Gp^Kiz) = PKz^ - cg{2pKz). (3.4) 

The calculation of the zeroes of Ip^K — equivalently, the global minimum points of Jp^xiz) — 
13 Kz^ — is greatly facilitated by the following observations made in Proposition 3.4 in [22]: 

1. The global minimum points of J(},k{z) — l3Kz^ coincide with the global minimum points 
of Gfi^K, which are much easier to calculate. 

2. The minimum values Yahizm{J i3,k{z) — f^Kz"^} and min^gK Gj3^K{z) coincide and both 
equal the canonical free energy K) defined in part (b) of Theorem 3.1. 

Item 1 gives the alternate characterization that 

^0,K = {-2 e [—1, 1] : z minimizes Gp^K{z)}- (3.5) 

In the context of Curie-Weiss-type models, the form of Gp^K is explained on page 2247 of [22]. 

As shown in the next two theorems, the structure of Efj^K depends on the relationship be- 
tween 13 and the critical value /3c = log 4. We first describe Ep^x for < /9 < /5c and then for 
13 > Pc- In the first case Sp^x undergoes a continuous bifurcation as K increases through the 
critical value Kc{(3) defined in (3.6); physically, this bifurcation corresponds to a second-order 
phase transition. The following theorem is proved in Theorem 3.6 in [22]. 

Theorem 3.2. For Q < (3 < (3c, we define 

For these values of (3, Sgx has the following structure. 

(a) ForO < K < K,{P), Sp^K = {0}. 

(b) For K > Kc{(3), there exists z{P, K) > such that Sp^x = {±z{(3, K)}. 

(c) z{(3,K) is a positive, increasing, continuous function for K > Kc{j3), and as K 
{Kc{pi))'^, z{P, K) — > 0. Therefore, Ep^x exhibits a continuous bifurcation at Kc{(3). 

For 13 e (0, /3c), the curve (/3, Kc{l3)) is the curve of second-order points. As we will see 
in a moment, for (3 e (/3c, 00) the BEG model also has a curve of first-order points, which we 
denote by the same notation (/3, Kc{l3)). In order to simplify the notation, we do not follow the 
convention in [22], where we distinguished between the second-order phase transition and the 
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first-order phase transition by writing Kc{(5) for < /3 < /3c as Ki \(5) and writing Kc{P) for 



We now describe Sg^K for f3 > (3c. In this case Sj3^K undergoes a discontinuous bifurcation as 
K increases tlirough an implicitly defined critical value. Physically, this bifurcation corresponds 
to a first-order phase transition. The following theorem is proved in Theorem 3.8 in [22]. 

Theorem 3.3. For all (3 > Pc, £(3,k htis the following structure in terms of the quantity Kc{(3), 
denoted by Kc^^ in [22] and defined implicitly for (5 > (5c on page 2231 of [22]. 



(a) ForO<K< Kc{(3), E^^k = {0}. 

(b) There exists z{(3, Kc{(3)) > such that Sp^K^ifS) = {0, ±z{l3, Kc{l3))}. 

(c) For K > Kc{(3) there exists z{(3, K) > Q such that E^^k = {^z{(3, K)}. 

(d) z{(3,K) is a positive, increasing, continuous function for K > KdP), and as K 
Kc{(3)^, z{P, K) —>■ z{(3, Kc{(3)) > 0. Therefore, E^^k exhibits a discontinuous bifurcation at 



We end this section by outlining the proofs of the laws of large numbers in (2.1) and (2.3) 
and its breakdown in (2.2). The upper large deviation bound in the LDP stated in part (a) of 
Theorem 3.1 implies that for any > and X > the limiting mass of Sn/n with respect 
to Pn,i3,K concentrates on the elements of E^^k. According to Theorems 3.2(a) and 3.3(a), 
E()^K '= {0} when < p < pc and < K '< K^p) and when p > pc and K < KdP). 
For these values of P and K, the laws of large numbers in (2.1) and (2.3) follow immediately. 
For P > Q and K > Kc{P), we have Eg^x = {±2;(/3, K)}, and so by symmetry the limit (2.2) 
follows. The proof of the limit (2.4) is postponed until after Theorem 4.2 because it requires 
more detailed information about the elements of E^^k when P > Pc and K — Kc{P). 

In the next section we present additional properties of the function Gp^x introduced in (3.4). 
These properties will be used in later sections to prove the scaling limits and the MDPs for 



p> PcasKl'\p). 



KciP). 



Sn/n 



1-7 



4 Properties of Gf3^K 



As we saw in (3.5), the global minimum points of 



Gp,k{z) = PKz^-cp{2PKz) 
= PKz'' -log\^^ 




1 + 2e-^ 



coincide with the elements of E/g^K, the set of equilibrium macrostates for the BEG model. In 
this section we study further properties of G/j,^ that will be used in later sections to prove the 
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scaling limits and the MDPs for Sn/ n}"^ with respect to Pn,f3,K and with respect to Pn,i3„,K„ for 
appropriate sequence (/?„, Kn) and for appropriate choices of 7. 

We first show that for any 7 G [0, 1) the P„^^ /^^-distribution of Sn/n^^'^ can be expressed 
in terms of Gp^^K„ and an independent normal random variable. The next lemma can be proved 
like Lemma 3.3 in [20], which applies to the Curie-Weiss model, or like Lemma 3.2 in [25], 
which applies to the Curie-Weiss-Potts model. In an equivalent form, the next lemma is well 
known in the literature as the Hubbard-Stratonovich transformation, where it is invoked to an- 
alyze models with quadratic Hamiltonians (see, e.g., [1, p. 2363]). After the statement of the 
lemma, we outline how we will use it in order to deduce the scaling limits of Sn/n^^"'. 

Lemma 4.1. Given a positive sequence {f3n,Kn), let Wn be a sequence of N{(}, (2/?„/C„)~^) 
random variables defined on a probability space (^2, Q). Then for any 7 G [0, 1) and any 
bounded, measurable function f 

f{x) ex.p[-nGf}^,Kn{x/n'^)]dx. 



f^expl-nGp^^K^x/n-y)] dx 

As we will see in Theorems 5.1, 6.1, and 7.1, the scaling limits have different forms de- 
pending on which of the following three sets K) lies in: the singleton set C containing the 
tricritical point (/3c, Kc{Pc))', the curve B of second-order points 

B = {{(5, K)eR':0<f3<f5c,K = K,{f3)}, 

and the single-phase region 

A = {((3,K) G < /3 < /3c,0 < X < K,((3)}. 

These sets are shown in Figure 1 in the introduction. 

We now indicate how we will use Lemma 4.1 to prove the scaling limits of Sn/n^^^ for 
7 G (0, 1/2]. Let {Pn, Kn) be a suitable positive sequence converging to (P , K) E A U B U C. 
Assume that Kn) and 7 are chosen so that the limit of the right hand side of (4.1) exists 
as n — > cxD. We first consider 7 < 1/2. Since /3„ and A'„ are bounded and uniformly positive 
over n, rewriting the limit of the left hand side in terms of characteristic functions shows that 
Wn/ n^/^"'"' does not contribute. Hence it follows that 

lim / f{SJn'-^)dPn,p^,Kr. (4.2) 

= lim — ; — ; — , ■ / f(x) exp[— k (x/rf)]dx. 

n-^^J^eM-nGp^,KA^/n^)]dxJJ^' ^"'^"^ ^ 



Costeniuc, Ellis, and Otto: Critical Behavior of Probabilistic Limit Theorems 



24 



From this formula we will be able to determine the scaling limits of Sn/ n} when Kn) — > 

{l3, K) e B Li C [Thms. 6.1, 7.1]. Using an analogous formula, we will be able to determine 
the MDPs of Sn/n'^-^ when (/5„, Kn) ^ {P,K) e BUC [Thms. 8.1, 8.3]. 

Now consider 7 = 1/2, which corresponds to the central-limit-type scaling for Sn in (2.5). 
In this case (4.1) yields 

lim / f{Sjn^^' + Wn)d{Pn,p„,K„xQ) (4.3) 

jRexp[-nG'^„,K„(a;/nV2)]cia; 

In contrast to when 7 G (0, 1/2), Wn now contributes to the limit. Again the use of characteristic 
functions enables one to determine the scaling limit of Sn/n}^'^ when Kn) K) e A 

[Thm. 5.1]. 

Formulas (4.2) and (4.3) suggest how to proceed in proving the scaling limits of Sn/n}"^. 
First consider (/?„, Kn) for which Gp^^K^ has a unique global minimum point at [Thms. 3.2(a), 
3.3(a)]. As (4.2) and (4.3) suggest, the forms of the scaling limits of Sn/n^^"' with respect to 
Pn..8n,K„ depend on the forms of the Taylor expansions of Gjs^^Kn in the neighborhood of the 
global minimum point 0. One of the attractive features of our analysis is that the same Taylor 
expansions can be used to handle sequences K^) for which Gp^^K^ has nonunique global 
minimum points. Such sequences arise naturally in the scaling limits and the MDPs to be proved 
later in the paper; in fact, it is precisely such sequences for which the MDPs yield the new class 
of distribution limits of the form (2.13) and (2.14). What makes it possible to use the same 
Taylor expansions regardless of the nature of the global minimum points of Gp^^K„ is Lemma 
4.4, the main technical innovation in this paper. 

Preliminary information on the forms of the relevant Taylor expansions is presented in Theo- 
rems 4.2 and 4.3. In the proofs of the scaling limits, in order to justify replacing nGp^^Knix/n'^) 
in (4.2) by n times the Taylor expansion evaluated atx/ri' and taking limits under the integral, 
one invokes the dominated convergence theorem, for which the appropriate bounding function 
depends on the particular sequence This will be handled on a case-by-case basis in 

subsequent sections. Finally, one must show that the contributions to the limit in (4.2) and (4.3) 
by all X for which x/n'^ lies in the complement of a neighborhood of is exponentially small. 
The relevant error estimate is given in part (c) of Lemma 4.4. Similar considerations apply to 
the proofs of the MDPs in section 8, for which the relevant error estimate is given in part (d) of 
Lemma 4.4. 

The steps outlined in the preceding paragraph for deducing the scaling limits of Sn/n^~'^ 
from (4.2) and (4.3) are well known in the related contexts of the Curie- Weiss model and the 
Curie-Weiss-Potts model. Scaling limits for these models are studied in [20, 21] and in [25] 
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for fixed values of the inverse temperature defining the corresponding canonical ensemble. In 

contrast to those earlier papers, our study of scaling limits for the BEG model necessitates a 
considerably more careful analysis because we work with the canonical ensemble Pn,/3„,i<:„> 
allowing sequences Kn) rather than only fixed values of {f3,K). 

The analysis of the Taylor expansions of in the neighborhood of a global minimum 
point involves the notion of the type of a global minimum point, which we next introduce. We 
temporarily consider any P > and any 7^ > and then specialize to {(3,K) E AU B U C. 
Let z be an element of £/3^k- Since G^^k is real analytic and 5 is a global minimum point, there 
exists a positive integer r = r{z) such that G^p^^{z) > and 

Gp,K{z) = Gp{z) + -^(^ - + 0{{z - ~zr+') asz^ ~z. 

We call r{z) the type of the global minimum point z. If r = 1, then = 2PK — 

{2(5KY{cp)"{2(5K~z), and if r > 2, then Gf\z) = -{2pKf^c^^'\z). 

In Theorem 6.3 in [22] the types of the elements of S^^k are determined for all /3 > 
and K > 0. In our study of scaling limits of Sn/'n}~'^ in the present paper, we focus on 
{l3,K) e AU B UC, for which S^^k = {0} [Thm. 3.2(a)]. Although the conclusion in 
[22] that for {(3,K) e B the type of equals 2 is correct, the formula for G^'^^(O) given in 
(6.6) in that paper has a small error. The correct formula for G^''^(0) is given in (4.10) with 

Theorem 4.2. For all {p,K) e AU BUG, 8,3,k = {0}. 

(a) For all {(3, K) E A, z = has type r = 1. 

(b) For all {p, K^{(3)) E B, z = has type r = 2. 

(c) For {(3, K) = l(3c, KdlSc)) EG,z = Ohas type r = 3. 

For all other values of P > and K > not considered in Theorem 4.2, the elements of 
£p^K all have type r = 1. This includes the values < P < Pc and K > Kc{P) [Thm. 3.2(b)] 
and the values P > Pc, K > [Thm. 3.3]. In these two cases the fact that the elements of Ep^K 
all have type r = 1 is proved in [22] in part (c) of Theorem 6.3 and in Theorem 6.4. 

We now point out how to prove the breakdown of the law of large numbers stated in (2.4), 
which holds for P > P^ and K = Kc{P). In this case, £j3^Kc{i3) = {0, ±2;(/9, X)}. Since 
each of the elements of Sp^KciP) ^^P^ r = 1, the limit in (2.4) is proved exactly as in part 
(c) of Theorem 2.3 in [25], which treats the breakdown of the law of large numbers for the 
Curie-Weiss-Potts model at /3 = Pc. In (2.4), 



(4.4) 
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where «o = [G'gJ,^(;3)(0)]-^/^ and «i = [G^^^^^fsMP, K^)]-'^'- 

We return to Lemma 4.1 and in particular to (4.2)-(4.3), which express the scaling limit of 
Sn/n^^^ in terms of the function nG i3^,kS'x / n^) ■ Using the information about the three differ- 
ent types of the global minimum point of Gjs^K at for K) e A, {j3, K) e B, and K) e 
C, we now indicate the three different forms of the Taylor expansion of nGp^^K^ {x/rO) needed 
to deduce the scaling limits of Sn/n^~'^. These involve the quantities G^^j^^(O), G^'^^^^(O), 

and G^^j^ (0), for the first two of which explicit formulas in terms of /3„ and Kn are given. 
As we will see in later sections, these formulas will guide us into how we should choose the 
sequences Kn) so that all the different scaling limits of Sn/n^^^ appear. Since G^„^k„ is 
symmetric around 0, all odd-order derivatives of this function evaluated at vanish; in addition, 

In order to state part (d) of the theorem, we define for > 

For < P < Pc this function coincides with the function Kc{P) defined in (3.6), while for 
/3 > (3c, K{(3) > KciP) [22, Thm. 3.8]. Thus for {/3, K) e B we have K = K^iP) = K{/3) 
while for {(3, K) e C we have /3 = /3c and X = Kc{pc) = K{pc). 

Theorem 4.3. Let Kn) be any positive bounded sequence and 7 any positive number. The 
following conclusions hold. 

(a) Assume that K„) {(3, K) e A. Then the type ofOE £p^K equals 1. In addi- 
tion, for any R > and for all x & M. satisfying \x/n'^\ < R there exists ^ — $^{x/n') e 
[—x/n', x/n'] such that 

nGp^,KS^/n^) = ^"'^^"^ x' + ^An{ax/n'^))x'. (4.6) 

The error terms An{Cix/n'^)) are uniformly bounded over n & N and x e {—Rn'^Rri'). 
Furthermore, as n ^ 00, Gg^^7^^(0) G^^^(O) > 0. 

(b) Assume that {j3n,Kn) {f3,Kc{(3)) G B. Then the type of E Sp^Kciis) '■^ 2. In 
addition, for any R > and for all a; e M satisfying \x/n'^\ < R there exists ^ = ^{x/n') e 
[—x/n'^x/n'] such that 



1 Gf J. (0) o 1 Gf J. (0) , 1 



The error terms Bn{^{x/7i')) are uniformly bounded over n G N and x G {—Rri^^RrC). 
Furthermore, as n ^ oc, Gg^;,„(0) ^ ^(^^.^(^^(O) = while G'^I^JO) ^ ^(^^.^(^^(O) > 0. 
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(c) Assume that {(3m — > {(3c, Kc{(3c))- Then the type ofO& Sp^^Kd/Sc) 3. In addi- 
tion, for any R > and for all x & M. satisfying \x/n''\ < R there exists ^ — ^{x/n') e 
[—x/v?, x/n'] such that 



nGp^^KS^/n^) = (4.8) 



The error terms Cn{i{x/n^)) are uniformly bounded over n G N and x G {—Rn',Rn'). 
Furthermore, as n ^ oo, Cjfj .^jO) ^ Gjfj,Ke(/3c) (0) = ^ <^?nUn(0) ^ '^PMPc^^) = ^ 

(d) For(3>0 we define K{(3) in (4.5). Then in (4.6)-(4.8) 

^(2) ... _ 2/3A(e^"+2-4/3„/r^) _ 2/3^/r4ir(/30 - K^] 

Gp.,KS^) - " K{^) ^^-^^ 



and 

G(V(0) = ^Pft;y(^ (4.10) 



Proof. In parts (a), (b), and (c) the type of the global minimum point at is specified in Theorem 
4.2. The formulas for G^^j_;^^(0) and G^^^ ;^^(0) in part (d) follow from an explicit calculation 
of the derivatives and from the formula for K{(3) given in (4.5). In addition, one evaluates the 
limits of the Taylor coefficients given in the last sentence of each part (a), (b), and (c) using the 
continuity of the derivatives G^^'j^(O) with respect to (3 and K and the fact that the type of the 
global minimum point of Gp^K at is, respectively, r = 1, r = 2, and r = 3. 

We now prove the form of the Taylor expansion given in part (c); the forms of the Taylor 
expansions given in parts (a) and (b) are proved similarly. By Taylor's Theorem, for any i? > 
and for all G IR satisfying |it| < R there exists ^ — ^{u)& [— xi, u\ such that 

^f}„,KnW- u H u H U +Cn[^[U))U, (4.11) 

where C„(^(it)) = G^^^ j^^{^{u))/7\. Because the sequence (/3„, Kn) is positive and bounded, 
there exists b G (0, oo) such that < Pn ^ b and < Kn < b for all n. As a continuous function 
of K, x) on the compact set [0, b] x [0, b] x [—R, R], Gg^j^{x) is uniformly bounded. It follows 

that G^p^ f^^{^{u)), and thus Cn{^{u)), are uniformly bounded over n G N and u G {—R, R). 
Multiplying both sides of (4.1 1) by n and substituting u = x/n? yields part (c). ■ 
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This completes our preliminary discussion of the Taylor expansions of nGp^^Kni^/n^) as 
they are needed to deduce the scaling limits of Sn/n^~^ via Lemma 4.1. In order to finalize our 
analysis of these scaling limits, we will have to prove that the contributions to the integrals in 
(4.2) and (4.3) by a; G M satisfying Ix/ri^l > R converge to as n — >^ cxd. In part (c) of the 
next lemma we prove that the convergence to is exponentially fast. The technical hypothesis 
in part (c) is satisfied in each of the theorems that proves the scaling limits [Thms. 5.1, 6.1, 
7. 1]. In part (d) of the next lemma we prove the exponentially fast convergence to of a related 
integral that arises in the proof of the MDPs. As we verify in the proof of Theorem 8.1, the 
technical hypothesis in part (d) is satisfied in that setting. The estimates in parts (c) and (d) 
are consequences of the LDP proved in part (b), which in turn follows from part (a) and the 
representation formula in Lemma 4.1. 

Lemma 4.4 is the main technical innovation in this paper. When adapted to the BEG model, 
the precursors of Lemma 4.4 given in Lemma 3.5 in [20] and Lemma 3.3 in [25] are able to 
handle only positive sequences K^) converging to {l3,K) e AU B U C for which Gp^^Kn 
has a unique global minimum point at 0. In order to handle sequences X„) for which 
Gi3„,K„ has nonunique global minimum points, the modifications that would be necessary in 
the precursors of Lemma 4.4 would introduce serious technical complications in the proofs of 
the scaling limits and the MDPs. By allowing us to handle any positive sequence [pn, Kn) 
converging io [13 , K) ^ A\J B \J C, parts (c) and (d) of Lemma 4.4 are universal bounds that 
enable us to avoid these technical complications altogether. 

Lemma 4.4. Let {j3n, Kn) be any positive sequence converging to {j3, K) G AVJ B VJ C and 
as in Lemma 4.1, let Wn be a sequence of N{0, (2/3„i^„)~^) random variables defined on a 
probability space {D,,T,Q). The following conclusions hold. 

(a) There exist Oi > and 02 > such that for all n & N and all x G R, Gp^^Knix) > 
ai{\x\ - ly - 02- 

(b) With respect to Pn,f3„,K„ x Q, Sn/n+ Wn/n}^'^ satisfies the LDP on M with exponential 
speed n and rate function Gp^x- 

(c) Given 7 > and R> 0,we define 



If the sequence yn is bounded, then there exists 03 > and 04 > such that for all sufficiently 



(d) Assume that there exist 7 > 0, i? > 0, m G (0, 1), 05 > 0, and Oe G M such that for all 




(4.12) 



large n 




exp[-nG^„,/<r„(a;/n^)] dx < a3exp(-na4) 0. 
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sufficiently large n 



Vr 



exp[-nG^„,/<:„(x/n^)] dx < 05 exp(n"a6). 



' {\x\<Rni} 

Then there exists 07 > such that for all sufficiently large n 



L 



exjg)[—nGjj^^Kn{^/n'^)\ dx < 2a^ex.]){—na^) — >• 0. 



{\x\>Rni} 

Proof, (a) Because the sequence (/J^j, Kn) is bounded and remains a positive distance from the 
origin and the coordinate axes, there exist < 61 < 62 < 00 such that 61 < /3n < ^2 and 
bi ^ Kn < 62 for all n e N. The conclusion of part (a) is a consequence of the elementary 
inequalities 

> PnKnx'' - 2p„Kn\x\ - log4 > bl{\x\ - If -hi- log 4. 
(b) We prove that for any bounded, continuous function ^ 



lim — log 

n— »oo n 



exp 



ml) 



n n 



1/2 



d{Pn,0^,K^ X g) = sup{V^(a;) - G0,k{.x)}. (4.13) 



This Laplace principle implies the LDP stated in part (b) [13, Thm. 1.2.3]. Gf^^K is continuous, 
and by part (a) of this lemma applied to the constant sequence Kn) — K), this function 

has compact level sets. Since {(5 , K) E AVJ B C , Gp^x has a unique global minimum point 
at 0, and therefore inf^gu ^^^/^(a;) = 0. It follows that Gp^x is a rate function. We now use 
Lemma 4. 1 with 7 = to rewrite the integral in the last display as 



exp 



nil) 



1 



d{Pn,fi^,K^ X Q) (4.14) 

/ exp[n{V'(x) - Gp^^K^x)}] dx. 
Jm. 



j^Qx^[-nGp^,KSx)]dx 

By part (a) of this lemma, there exist M > and ag > having the following three 
properties: 



1- Gji^^K^i^) > asx^ for all n e N and all x e R satisfying \x\ > M. 



2. The supremum of ^ — Gj3^K on M is attained on the interval [— M, M]. 
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3. Let A = sup^^^{7P{x) - Gp,k{x)}. Then \\^\\^ - < -|A| - 1. 

Since Gp^^Kn converges uniformly to Gp^K on [— M, M], we have for any 5 > and all suffi- 
ciently large n 



exp(—nS) / exp[n{'0(x) — Gp^xix)}] dx 

J{\x{<M} 

< ex.p[n{ilj{x) -Gf)^^KAx)}]dx 

J{\x\<M} 

< exp(n5) / exp[n{'ilj{x) — Gp^K{x)}] dx. 

JUx\<M\ 



In addition, by items 1 and 3 

/ exp[n{^(a;) - G/3„,i^„(a;)}] dx 

J{\x\>M} 

< exp[n||'0||oo] / exp[—nasx^] dx 

Jux\>M\ 



< -4r- cxp[n||V'||oo - nasM^] 

nlvl as 



We now put these estimates together. For all sufficiently large n we have 

exp(—nS) / exp[n{t/j{x) — Gp^xi^)}] dx 

J{\x\<M} 

< / exp[n{ijj{x) - Gp„,K„{x)}]dx 

Jm. 

< exp(nS) [ exp[n{ijj(x) — G3 K{x)}]dx -\ exp[— n(|A| + 1)]. 

J{\x\<M} ' riMas 

Since by item 2 

lim —log / exi)[n{^{x) — Gp xix)}] dx 



'{\x\<M} 

sup {i>{x) - Gfj^xix)} = sup{V^(x) - Gp^xix)}, 

{12;|<M} xeR 
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we see that 



< liminf - log / cxp[n{ijj{x) - G/3„,x„(a;)}] dx 

n-*oo n 

< limsup - log / exp[n{'0(x) - Gp^^Kn{x)}] dx 

< sup{ilj{x) - Gp,k{x)} + S, 
xeM. 

and since 5 > is arbitrary, it follows that 

lim - log / ex.p[n{ijj{x) - Gp^^Kn{x)}] dx = sup{ip{x) - Gp^xix)}. 

n-»oo n x&t 

We combine this limit with the same limit for = and use (4.14) together with the fact that 

ini^^^Gp^Kix) = G'/3,A-(0) = 0, concluding that 



lim — log / exp 



n— >c» n 



A"xf2 



— + ^2 



= sup{^(a;) - Gp^k{x)} - inf Gp^xix) = sup{^(x) - Gf^^Kix)}. 

This is the Laplace principle (4.13). The proof of part (b) is complete. 

(c) Since G^^k has a unique global minimum point at 0, the LDP proved in part (b) implies 
the existence of ag > such that for all n e N 



Pn,(^r.,Kr. X + ^ ^ (-^,^)} < exp(-na9). 

Using Lemma 4.1, we rewrite the probability in the last display as 

Pn,,„,K„Xg|^ + ^^(-i?,i?)} 

exp[-nG;3„,if„(x/n^)] dx 



(4.15) 



(4.16) 



J{\x\>Rn-'} 
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where is defined in (4.12) and 

Zn= exp[-nGf3^,K„{x/n'^)]dx. 

J{\x\>Rn')} 

Since by hypothesis the sequence is bounded, there exists y > such that yn < y for all n. 
It follows from (4.15) and (4.16) that for all sufficiently large n 

\zn < Zn{l - exp(-na9)) < |/„exp(-ria9) < i/exp(-na9) 

and thus for all sufficiently large n, Zn < 2yexp(—nag). This completes the proof of part (c). 
(d) Exactly as in the proof of part (c), we have for all sufficiently large n 

\z„ < Znil - cxp{-nag)) < i/„cxp(-na9). 

Since by hypothesis y„ < exp(n"a6) and u e (0, 1), it follows that for all sufficiently large n 

Zn < 2a5exp(— nog + n^oe) < 2a5 exp(— ^09/2). 

This completes the proof of part (d). ■ 

In the next section we begin our analysis of the scaling limits of Sn/n^~^ in the simplest 
case by considering (/?„, Kn) {(3, K) E A. In the two sections following the next one, we 
will uncover a wider variety of scaling limits by considering sequences K^) converging to 

iP, K,{i3)) e S and to {13,, K,{P,)) e C. 

5 1 Scaling Limit for Kn) {l3, K) e A 

In this short section, we deduce the unique scaling limit of Sn/n^~"' when {/3n, Kn) is any 
positive sequence converging to (/?, K) E A. The unique global minimum point of G^^k at 
has type r = 1 [Thm. 4.2(a)]. As the next theorem shows, the scaling limit with respect 
to Pn,/3„,Kn has the form of a central limit-type theorem that is independent of the particular 
sequence chosen. In addition, the only value of 7 for which Sn/ n^~'^ has a nontrivial limit is 
7 = 1/2. We are including this scaling limit in order to highlight the much more complicated 
behavior of the scaling limits of Sn/n^~'^ in the subsequent two sections, in which Kn) — > 
(/?, Kc{(3)) E B and Kn) (/5c, K,{l3c)) E C and in which different forms of the limit 
can be obtained by choosing different sequences. 

The following theorem, stated for < /5 < /5c and < K < Kc{P), is also valid for (3 > (3, 
and Q < K < Kc{(3), and the proof is essentially the same. The key observation is that for 
13 > /5c, we have K{(3) = {e^ + 2)/(4/5) > K,{(3) [22, Thm. 3.8]. Hence if K < K,{p), then 
also K < K{(3) and thus G^^lO) in (5.2) is positive. 
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Theorem 5.1. Let Kn) be an arbitrary positive sequence that converges to {(3, K) e A; 
thus (5 and K satisfy < (3 < (3c cind Q < K < Kc{(3). Then 

Pn,ii.a,Kr,{Sn/n^^'^ e dx} =^ exp(-C2X^) dx, 

where C2 > is defined by 

C2 = J • = (3[K{(3) - K]. (5.1) 

Thus the limit is independent of the particular sequence {Pn, -^n) ^hat is chosen. 

Proof. We use the Taylor expansion in part (a) of Theorem 4.3 with 7 = 1/2. By continuity, 
G^^j (0) given in (4.9) converges to 

which is positive since < K < Kc{(3) = K{(3). For any i? > the error terms An{x/n^^'^) in 
the Taylor expansion are uniformly bounded over n G Nanda; e {-Rn^/^ , Rn^/"^). It follows 
that for all X e R 

lim nG0^^KSx/n"')^\G%{Q)x' 

and that R > can be chosen to be sufficiently small so that for all sufficiently large n and all 
X e R satisfying \x/n^^'^\ < R 

Since exp[— G^^^(0)a;^/4](ia; < 00, the dominated convergence theorem implies that for any 
bounded, continuous function / 

lim / fix) e^p[-nGp^,KMn'/')]dx^ [ /(x) exp[-G'gi,(0)xV2] dx. 

The existence of this limit implies that the sequence i/n = /r|a;|<ij„i/2| cxp [—nGp^^Xn {x / 'n}^'^)]dx 
is bounded. Hence, combining this limit with part (c) of Lemma 4.4 yields 

lim / f{x)eM-nGf,„,KMn'/')]dx= [ /(a:) exp[-Ggi,(0)a:V2] c^x. 



Costeniuc, Ellis, and Otto: Critical Behavior of Probabilistic Limit Theorems 



34 



We now augment this limit with the same Umit for / = 1 and use (4.3) to obtain 

lim / f{SJn'/' + W^)d{P^,f,,K.^f,)XQ) 

We omit the straightforward argument using characteristic functions that enables one to deduce 
from the last display that 

^n,/3„,K„{5'„/n^/^ e dx} =^ exp{-C2X^) dx, 

where C2 is given by the first equality in (5.1). A similar argument involving moment generating 
functions is given on pages 70-71 of [25]. The positivity of C2 and the second formula for C2 
given in (5.1) follow from (5.2). This completes the proof of the theorem. ■ 

In Theorem 8.2 we prove an MDP for Sn/rt^^'^ that is related to the scaling limit proved in 
Theorem 5.1. As in the latter theorem, the form of the MDP is independent of the particular 
sequence (/?„, Kn) converging to K) E A. In the next section we see the first example 
of scaling limits for Sn/n^'"' where different forms of the limit can be obtained by choosing 
different sequences Kn) — > (/?, Kc{P)) G B. 



6 4 Scaling Limits for K^) Kc{(3)) G B 

In this section we determine the scaling limits of Sn/n^~"' with respect to Pn,8„,K„, where 
(/?„, Kn) is an appropriate positive sequence converging to (/?, /(',.(/?)) G B and 7 G (0, 1/2) is 
appropriately chosen. We recall that B is the curve of second-order points for the BEG model. 
For any K) G B, we have < /? < /?c = log 4 and 

The scaling limits that we obtain involve limiting densities proportional to exp[— G(a;)], 
where G takes one of the 4 forms of an even polynomial of degree 4 or 2 satisfying G{0) = 
and G{x) — > 00 as \x\ — > 00. There are 3 such G's of degree 4; namely, G{x) = c^x'^, where 
C4 > and G{x) — kpx'^ + C4X^, where C4 > and either A; > or A; < 0. There is also 1 such 
G of degree 2; namely, G{x) — kpx'^, where A; > 0. These 4 cases are all obtained in Theorem 
6. 1 ; the forms of the limits depend on the choice of K^ — > Kc{P) but are independent of the 
choice of (3n ^ 
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In order to determine the forms of the scaling limits of Sn/n^~^ with respect to Pn,/3„,i<:„, 
we start by recalling the Taylor expansion given in part (b) of Theorem 4.3. For any 7 > and 
i? > and for all x e M satisfying \x/rf\ < R there exists ^ e [—x/rf, x/rf] such that 



1 G^^^r, (0) , 1 G^^^r, (0) , 1 



The error terms Bn{i(x/n')) are uniformly bounded over n G N and x G {—Rrf^Rrf). 
According to part (b) of Theorem 4.2, the unique global minimum point of Gp^KciP) ^ 
type 2. Hence by continuity, as n — > 00, 



while G^pIk^{Q) G'^\(^)(0) > 0. We recall that in the last display K{I3) = (e'^ + 2)/(4/3) 
for /? > 0"' " 

Fixing /3 G (0, j3c), we let /9„ be an arbitrary positive sequence that converges to (3, and we 
let ^ be a positive number. The key insight is to choose so that G^^^ (0) ~^ at a rate 
where counterbalances the term l/n^'^~^ appearing in (6.1). Since 2PnKn/K{Pn) 
has the positive limit 2/3 as n —> 00, we achieve this by choosing k and defining 

Kn = K{l3n) - kin'. (6.2) 

Since Pn ^ P and K{-) is continuous, it follows that Kn — > K{P) — Kc{P). Hence 

gZkM = ■ ^ = ^ ■ C^'. »">e„ ) > and ^ 2f,. 

With these choices (6.1) becomes 

1 kr^'^^ 1 G^^^ fo) 1 

-G,.,K„(a;/n-) = + ^^^^^-^ + (6.3) 

As we will see in Theorem 6.1, the scaling limits depend on the value of 7 and on through 

the value of 9, but are independent of the sequence Pn ^ P- 

In the last display we assume that the coefficients multiplying and both appear with 
nonnegative powers of n and that at least one of these two coefficients has n to the power 0. 
Then in the limit n — > cxd any coefficient including the error term that has a positive power of n 
will vanish while any coefficient that has n to the power will converge to a positive constant. 
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This preliminary analysis shows the possibility of multiple scaling limits for different choices 
of 7 and 9. In order to confirm this possibility, we define 

V = min{27 + 6* - 1, 47 - 1} 

and focus on the cases in which i; = 0. As we will see in the final section of the paper, v < 
corresponds to 4 different MDPs for Sn/n^^^- On the other hand, if f > 0, then one obtains 
neither scaling limits nor MDPs. 

In the next theorem we show that v = corresponds to 3 different choices of 7 and 9, which 
in turn correspond to 4 different sequences Kn in (6.2). The additional sequence arises because 
when is not the highest order term in the scaling limit (cases 3-4), k can be chosen to be 
either positive or negative. As shown in Table 6.1 in part (b) of the theorem, for each of these 
4 different sequences we obtain 4 different scaling limits of Sn/n^~'^. In case 1 we can also 
choose k to be any real number; this affects only the definition of the sequence Kn, not the form 
of the scaling limit. 




7 

1/4 1/2 

Figure 4: Influence of A and B on scaling Limits when (/?„, Kn) (/?, KdP)) G B 



The results of the theorem confirm one's intuition concerning the influence of the regions 
on the scaling limits. Of the 4 cases, case 1 corresponds to the largest values of — namely, 
9 > 1/2 — and thus the most rapid convergence of Kn — ^ Kc{P). In this case only B influences 
the form of the limiting density, which is proportional to exp(— C4X^); C4 defined in (6.5) is 
positive since < e^" = 4. By contrast, case 2 corresponds to the smallest values of — 
namely, 9 E (0, 1/2) — and thus the slowest convergence of Kn Kc{P). In this case only A 
influences the form of the limiting density, which is proportional to exp(— /9x^); thus we have 
Sn/n^''^ converging in distribution to a normal random variable even though the non-classical 
scaling is given by n^~'^ , where 7 = (1 — 6')/2 G (1/4, 1/2). Finally, cases 3 and 4 correspond 
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to the critical speed 9 — 1/2. In this case both A and B influence the form of the Umiting 

density, which is proportional to cxp(— A';/?x^ — c^x^^) with C4 > and either A: > or A; < 0. In 
Figure 4 we indicate the subsets of the positive quadrant of the 9-^ plane leading to the 4 cases 
just discussed. Using Table 5.1, one easily checks that as 9 increases through the critical value 
1/2, the values of 7 in the scaling limit change continuously while the forms of the limiting 
densities change discontinuously. 

Theorem 6.1. For fixed /? e (0, /3c), let [3n be an arbitrary positive sequence that converges to 
(5. Given 9 > and k ^ 0, define 

= K{I3^) - kln\ 

where K{(3) = (e'' + 2)/{A(3)for (3 > 0. Then X^) ^ {(3, K^{(3)) e B. Given 7 e (0, 1), 
we also define 

G{x) = S{v, 2j + 9- l)kj3x^ + 5{v, 47 - l)c4a:^ (6.4) 
where d{a, b) equals lifa = b and equals Oifa^b and C4 > is given by 

^ _ ggk(/3)(0) _ mKm\^ - e') _ (e^ + 2)^(4 - e^) 
' 4! 4!(e^ + 2)2 23-4! ' 

The following conclusions hold. 

(a) Assume that v — min{27 + ^ — 1,47 — 1} equals 0. Then 

Pn,pn,Kn{Sn/n^~^ e dx} =^ exp[-G(a;)] dx. (6.6) 

(b) We have v = if and only if one of the 4 cases enumerated in Table 6.1 holds. Each of 
the 4 cases corresponds to a set of values of 9 and 7, to the influence of one or more sets B and 
A, and to a particular scaling limit in (6.6). In case 1 the choice ofk^M. does not affect the 
form of the scaling limit. 



case 


values of 9 


values of 7 


scaling limit of S^/n^ 


influence 








1 


9>i 


7 = i 


exp(— 042;^) dx 


B 






C4 > 0, A; e R 


2 




7 ~ ~2~ ^ (4' 2) 


exp(— dx 


A 






k>0 


3-4 




T=i 


exp(— A;/3a;^ — c^x^) dx 


A + B 






k^O 



Table 6. !.• Values of 9 and 7 and scaling limits in part (b) of Theorem 6.1 
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Note. Let Pn — P for all n. The constant sequence Kn) — {P, Kc{(3)) for all n corresponds 
to the choice 9 = oom case 1. As in the proof of case 1, one shows that Pn,0,Kc{()){Sn/'n^~^^^ £ 
dx) =^ ex.^{—CiX^)dx. This scaling limit was mentioned in (2.6). 

Proof of Theorem 6.1. We first prove part (b) assuming part (a), and then we prove part (a), 
(b) V — min{27+^ — 1, 47— 1} equals if and only if each of the quantities in this minimum 

is nonnegative and one or more of the quantities equals 0. As (6.4) makes clear, 47 — 1 = 
corresponds to the influence of B and 27 + — 1 = to the influence of A. We have the 
following 4 mutually exclusive and exhaustive cases, which correspond to the 4 cases in Table 
6.1. 

• Case 1: Influence of B alone. 27 + — 1 > 0, 47 — 1 = 0, and /c e M. In this case 
7 = 1/4 and ^>1 — 27=1/2, which corresponds to the second and third columns for 
case 1 in Table 6.1. 

• Case 2: Influence of A alone. 27+^-1 = 0,47-1 > 0, and/c > 0. In this case 7 > 1/4 
and ^ = 1 - 27 < 1/2. Since 9 must be positive, we have 7 = (1 - e)/2 e (1/4, 1/2). 
Hence case 2 corresponds to the second and third columns for case 2 in Table 6.1. 

• Cases 3-4: Influence of A and B. 27 + 6* — 1 = 0, 47 — 1 = 0, /c > for case 3, and 
A; < for case 4. In these 2 cases 7 = 1/4 and ^^ = 1 — 27 = 1/2, which corresponds to 
the second and third columns for cases 3 and 4 in Table 6.1. 

In cases 1, 2, 3, and 4 we have, respectively, G{x) = c^x^, G{x) = k/3x'^ with k > 0, G{x) = 
kfix"^ + c^x^ with A; > 0, and G{x) = kdx^ + C/^x'^ with A: < 0. In combination with part (a), 
we obtain the 4 forms of the scaling limits listed in the last column of Table 6.1. 

(b) We prove the 4 scaling limits corresponding to the 4 cases listed in Table 6.1. As the 
discussion prior to the statement of the theorem indicates, the quantity v — min{27+^ — 1, 47 — 
1} is defined in such a way that in each of the 4 cases defined by the choices of 9, 7, and k in 
Table 6.1, we have for each x e M 

lim nGp^^KS^/n') = G{x). 

n—foo 

Since in each case we have 7 e [1/4, 1/2), the term Wn/n}-^'^^^ in (4.1) does not contribute to 
the limit n ^ 00. Hence we can determine the scaling limits of Sn/n}"'^ by using (4.2). In 
order to justify taking the limit inside the integrals on the right hand side of (4.2), we return to 
(6.3) and use the fact that for all sufficiently large n, gP > and Gjfj,^„(0) > 0. It follows 
that R > can be chosen to be sufficiently small so that for all sufficiently large n and all x e M 
satisfying \x/n'^\ < R there exists a polynomial H{x) satisfying 

nGf,^,KMn'^)>H{x) (6.7) 
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and Jjg exp[—H{x)]dx < oo. In case 1 when /c > as well as in cases 2 and 3, H{x) — G{x) /2; 
in case 1 when A; < and in case 4, which corresponds to k < 0, 

H{x) = -2\k\l3x'^ + 04x72. (6.8) 

The last two displays in combination with the dominated convergence theorem imply that for 
any bounded, continuous function / 



lim / 



{\x\<Rni} 



f{x)exp[-nGp„,Kr,{x/n'^)]dx= / f{x)exp[-G{x)]dx. 

Jr 



The existence of this limit implies that the sequence yn = J^IxIkRui} Pn,Kr,{x /n^)]dx 

is bounded. Hence, combining this limit with part (c) of Lemma 4.4 yields 

lim / f{x) exp[-nGi3^^Kn{x/n'^)]dx = / f{x)ex.p[-G{x)]dx. 
Jr Jr 

If we augment this limit with the same limit for / = 1 and use (4.2), then we conclude that in 
each of the 4 cases 

lim / f{Sjn'~^) dPr,,p^,K. = ^ ■ ! fix) exp[-G(x)] dx. 

n^oojj^„ I exp[-G{x)\dx 7k 



J^exp[-G{x)]dx 

This yields the scaling limits in part (a). The proof of the theorem is complete. ■ 

This finishes our analysis of scaling limits for Sn/n^~"' with respect to Pn,f3ri,Kn' where the 
sequence Kn) converging to Kc{P)) G S is defined in Theorem 6.1. This analysis is 
a warm-up for the even more interesting analysis of the scaling limits for sequences X^) 
converging to the tricritical point. 



7 13 Scaling Limits for {f3n, K^) -> (/3e, K^{l3c)) 

In Theorem 6.1 we obtained 4 forms of scaling limits for Sn/n^~"' using sequences Kn) 
converging to a second-order point {(3, Kc{(3)) G B. The limiting densities are proportional to 
exp[—G{x)], where G takes of the 4 forms of an even polynomial of degree 4 or 2 satisfying 
G{0) — and G{x) — > oo as — > oo. In each case the form of the limit is independent of 
the choice of /9„ (3 but depends on the choice of Kn — > Kc{P). Like the BEG model at 
{P, Kc{P)) E B, the Curie-Weiss model has a second-order phase transition at a critical inverse 
temperature f3c. The 4 scaling limits and the 4 MDPs analyzed in Theorem 8. 1 are analogous to 
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the scaling limits and MDPs that hold in the Curie- Weiss model when the inverse temperature 

converges to Pc along appropriate sequences /9„ [14]. However, the 13 scaling limits proved in 
the present section and the 13 analogous MDPs obtained in Theorem 8.3 depend on the nature 
of the tricritical point, a feature not shared with the Curie-Weiss model. 

We now use the insights gained in the preceding section to study the more complicated 
problem of scaling limits for Sn/n^~^ using sequences K^^ converging to the tricritical 
point (/3c, Kc{(3c)) — (log 4, 3/[2 log 4]). As in the preceding section, we choose ^ > 0, /c 7^ 0, 
and 

Kn = K{Pn)-k/n', (7.1) 

where K{I3) = [e^ + 2)/ (4/5) for /3 > 0. In contrast to the preceding section, we now also have 
to pick the sequence /5„ appropriately. Theorem 7.1 shows that 13 scaling limits arise for differ- 
ent choices of 9, 7, and the parameter appearing in the definition of /?„. The limiting densities 
are proportional to exp[— G(a;)], where G takes one of the 13 forms of an even polynomial of 
degree 6, 4, or 2 satisfying G'(O) = and G{x) ^ 00 as \x\ 00. 

In order to determine the forms of the scaling limits for Sn/n}~'^ with respect to Pn,f3„,Kn'> 
we use the Taylor expansion given in part (c) of Theorem 4.3. For any 7 > and R > and 
for all X e ]R satisfying Ix/ii^l < R there exists ^ e [—x/rf, x/ii^] such that 



nGp^,KA^/n^) = (7.2) 
1 G^?j,(Q) , 1 Gfj,(0) , 1 Gff,(0) , 1 

„27-l 9 ^47-1 A\ ^67-1 R ^77-1 "^^^ ' " 



The error terms C„(^(a;/n^)) are uniformly bounded over n e N and x G {—Rrf,Rn^). 
According to part (c) of Theorem 4.2, the unique global minimum point of G/3^^Kc(i3c) ^ has 
type 3. Hence by continuity, as n — > 00, 



.(2) _ 2(3^K^[K(j3rl- K^] ^(2) 

K{l3n) 



GfK (0) = ^ G^^)^ (0) = 0, 



while<^„(0)-G(^)^^(,^)(0) = 2.3^ 

As in the preceding section, we choose Kn as in (7.1) so that G^^^ (0) ^ at a rate 
where counterbalances the term l/n^^'^^ appearing in (7.2). We also choose /3n so that 
Gf^ Kn(^) — >^ at a rate 1/n", where l/n°' counterbalances the term l/n^'^~^ appearing in 
(7.2). This is achieved by choosing a > and either 6 > or 6 < and then defining /3„ by the 
logarithmic formula 

/3n = log(4 - = log(e^^ - 6/^"); (V.3) 
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if 6 > 0, then /3„ is well defined for all sufficiently large n. Since and K{-) is 

continuous, it follows that Kn) {(3c, Kc{(3c)). With this choice of Kn) we have 

g'ZkM = -,-^M^ = -,-C!i}\ where C^) ^ 2/?„ (7.4) 



and 



Gf, (0) = ± . f/"^:> = ± . CW. „he„ ^ 4 > 0. (7.5) 

The dependence of Kn) in (7.1) and (7.3) upon a and ^ is complicated; because is a 
function of a, is both a function of and, through /9„, a function of a. However, the a and 
9 decouple nicely when (7.4) and (7.5) are substituted into (7.2), yielding 



nGp^,KS^/n^) (7.6) 
- ^2^+6-1 2! + n47+a-i 4! + 6! + ^^'^-^^W^ • 



We continue the analysis as in the preceding section. Let us suppose that in the last display 
the coefficients multiplying x^, x*, and all appear with nonnegative powers of n and that at 
least one of the coefficients has n to the power 0. Then in the limit n ^ oo any coefficient 
including the error term that has a positive power of n will vanish while any coefficient that has 
n to the power will converge to positive constants. In order to analyze the various cases, we 
define 

'u; = min{27 + 61 -1,47 + a -1,67-1}, (7.7) 

and focus on the cases in which w = 0. As we will see in the final section of the paper, w < 
corresponds to 13 different MDPs for Sn/n}''^. On the other hand, if w > 0, then one obtains 
neither scaling limits nor MDPs. 

In the next theorem we show that w = corresponds to 7 different choices of 7, 9, and 
a, which in turn correspond to 13 different sequences (/5„, Kn) defined in (7.1) and (7.3). The 
additional sequences arise because when x'^ is not the highest order term in the scaling limit 
(cases 4-5, 8-13), b can be chosen to be either positive or negative; similarly, when x"^ is not 
the highest order term in the scaling limit (cases 6-13), k can be chosen to be either positive 
or negative. As shown in Table 7.1 in part (b) of the theorem, for each of these 13 different 
sequences we obtain a different scaling limit of Sn/'n}~'^. 

The limiting densities in cases 1, 4-7, and 10-13 are new. In cases 2, 3b, 8, and 9 we obtain 
the same forms of the limiting densities as in Theorem 6.1, where we considered Kn) — > 
(/?, K) e B. However, the values of 7 in the corresponding scaling limits in the two theorems 
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are different. By contrast, the values of 7 and 9 as well as the forms of the limiting densities are 
the same in case 3a in Theorem 7.1 and in case 2 in Theorem 6.1. 

There are yet further possibilities concerning the sign of b and k. In all the cases in which 
no x"^ term appears in the scaling limit (cases 1, 3, 6, 7), we can choose b to be any real number. 
Similarly, in all the cases in which no x'^ term appears in the scaling limit (cases 1, 2, 4, 5), we 
can choose k to be any real number. Although the choice of ft or A; affects the definition of the 
sequence Kn), it does not affect the form of the scaling limit. 

Through the terms x^, x^, and appearing in the limiting densities, the scaling limits 
correspond to the influence of one or more of the sets C, B, and A. The influence of the 
various sets upon the form of the scaling limits is shown in Figure 2 in the introduction, and 
details are given in Table 7.1, which is included in part (b) of the next theorem. Case 3, which 
corresponds to the influence of A alone, has two subcases, labeled 3a and 3b in Table 7.1. Case 
3 a corresponds to the lower region labeled A in Figure 2 and case 3b to the upper region labeled 
A in Figure 2. Using Table 7.1, one easily checks that as (a, 9) crosses any of the lines in Figure 
2 labeled A + B, A + C, or B + C, the values of 7 in the scaling limits change continuously 
while the forms of the limiting densities change discontinuously. 

Theorem 7.1. Given a > 0, 9 > 0, b ^ 0, and k ^ 0, define 

Pn = log(4 - b/n'') = log(e''= - b/n'') and = K{Pn) - k/n\ 

where K{p) = (e^ + 2)/(4/?)/or l3 > 0. Then K^) (/5c, i^c(/3c)). Given 7 e (0, 1), we 
also define 

G{x) = S{w, 27 + ^ - l)kpcx'^ + d{w, 47 + a - l)6c4X^ + d{w, 67 - l)c6X^, (7.8) 

where C4 = 3/16 and cq — 9/40. The following conclusions hold. 

(a) Assume that w — min{27 + ^ — 1,47 + 0; — 1,67 — 1} equals 0. Then 

Pn,i3^,KASn/n'-^ G dx} =^ exp[-G(a;)] dx. (7.9) 

(b) We have w = if and only if one of the 13 cases enumerated in Table 7.1 holds. Each of 
the 1 3 cases corresponds to a set of values of 9, a, and 7, to the influence of one or more sets 
C, B, A, and to a particular scaling limit in (7.9). The form of the scaling limit is not affected 
by the choice ofb&M. in cases 1, 3, 6, and 1 and by the choice ofk&M. in cases 1, 2, 4, and 5. 
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case 


values of a 


values of 7 


scaling limit of 5„/n^ 


influence 


values of 9 


1 


a > i 




exp(— cgx^) dx 
c6>0,beR,keR 


C 


^ > 1 


2 


ae (0,i) 




exp(— 6C4X'') dx 

C4 > 0, 6 > 0, e R 


B 


e>^ 


3a 




7 ~ ^ (4) 2) 


exp{—kPcx'^) dx 
k>0,beR 


A 


9 e (0, A) 


3b 






exp{—k/3cx'^) dx 
k>0,beR 


A 


a > 29-1 


4-5 






exp(— 6c4X^ — cqx^) dx 
b^O,keR 


B + C 




6-7 


a > \ 




exp(— /c/JcO;^ — CqX^) dx 
kj^O,beR 


A + C 





8-9 


ae (0,i) 




exp(— /c/3cX^ — bc4x'^) dx 
ky^O,b>0 


A + B 


^ 2 ^ V2' 3/ 


10-13 


a = 1 


7=i 


exp{—k(3cx'^ — 6c4X^ — cqx^) dx 
A; ^ 0, 6 


A + B + C 






Table 7.1; Values of a, ^, and 7 and scaling limits in part (b) of Theorem 7.1 



Note. The constant sequence {Pn,Kn) = {Pc,Kc{Pc)) for all n corresponds to the choices 
a = ^ = 00 in case 1. As in the proof of case 1, one shows that Pn,/3c,Kc{f3c){Sn/n^~^^^ G 
dx} =^ exp{—CQX^)dx. This scaling limit was mentioned in (2.7). 

Proof of Theorem 7.1. We first prove part (b) from part (a) and then prove part (a). 

(b) w = min{27 + ^ — 1, 47 + a — 1, 67 — 1} equals if and only if each of the quantities 
in this minimum is nonnegative and one or more of the quantities equals 0. As (7.8) makes 
clear, 67 — 1 = corresponds to the influence of C, 47 + « — 1 = to the influence of 
B, and 27 + 6* — 1 = to the influence of A. We have the following 13 mutually exclusive 
and exhaustive cases, which correspond to the 13 cases in Table 7.1. In each of the cases the 
equalities and inequalities expressing the influence of one or more sets C, B, and A are easily 
verified to be equivalent to the equalities and inequalities involving a, 9, and 7 given in the 
second and third columns of Table 7.1. Case 3, the most complicated, divides into two subcases 
depending on the value of a. 



• Case 1: Influence of C alone. 27 + 61 - 1 > 0, 47 + a - 1 > 0, 67 - 1 = 0, 6 e R, and 
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keR. 

• Case 2: Influence of B alone. 27 + 6' - 1 > 0, 47 + a - 1 0, 67 - 1 > 0, 6 > 0, and 

keR. 

• Case 3: Influence of A alone. 27 + 6^ - 1 = 0, 47 + a - 1 > 0, 67 - 1 > 0, A; > 0, and 

beR. 

• Cases 4-5: Influence of B and C. 27 + ^ - 1 > 0, 47 + a - 1 = 0, 67 - 1 = 0, 6 > 

for case 4 and 6 < for case 5, and k eR. 

• Cases 6-7: Influence of A and C. 27 + ^ - 1 = 0, 47 + a - 1 > 0, 67 - 1 = 0, A; > 

for case 6 and A; < for case 7, and 6 e M. 

• Cases 8-9: Influence of A and 27 + ^ - 1 = 0, 47 + a - 1 = 0, 67 - 1 > 0, A; > 
for case 8, /c < for case 9, and b > 0. 

• Cases 10-13: Influence of A , and C 27 + ^ - 1 = 0, 47 + a - 1 = 0, 67 - 1 = 0, 

k > and 6 > for case 10, A; < and 6 > for case 11, A; > and 6 < for case 12, 
and A; < and 6 < for case 13. 

In each of the 13 cases the form of G{x) follows from (7.8). In combination with part (a), 
we obtain the 13 forms of the scaling limits listed in the last column of Table 7. 1 . 

(a) The proof of the 13 scaling limits follows precisely the pattern of the proof of the 4 
scaling limits listed in part (b) of Theorem 6.1. As the discussion preceding the statement of 
Theorem 7.1 indicates, the quantity w — min{27 + ^ — 1, 47 + a — 1, 67 1} is defined in 
such a way that in each of the 13 cases defined by the choices of a, 9, 7, k, and b in Table 7.1, 
we have for each x e R 

lim nGp^^Knix/n^) = G{x). 

n— »oo 

Since in each case we have 7 e [1/6, 1/2), the term W^/n^/^"'"' in (4.1) does not contribute to 
the limit n — > 00. Hence we can determine the scaling limits of Sn/n^~'^ by using (4.2). In 
order to justify taking the limit inside the integrals on the right hand side of (4.2), we return to 
(7.6) and use the fact that for all sufficiently large n, C^^ > 0, Ck^^ > 0, and G^^^^^JO) > 0. It 
follows that R > can be chosen to be sufficiently small so that for all sufficiently large n and 
all X e R satisfying Ix/rf] < R there exists a polynomial H(x) satisfying 

nG0^,KMn^)>H{x) (7.10) 

and cxp[— if (a;)] < cxo. We define H{x) — G{x)/2 in all the cases in which both b > and 
A; > (cases 1-4, 6, 8, 10). Otherwise, a suitable polynomial H can be found as in (6.8); the 
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details are omitted. As in the proof of Theorem 6.1, the dominated convergence theorem and 
part (c) of Lemma 4.4 imply that for any bounded, continuous function / 

lim / f{x)exp[-nGfs„,K„{x/n'^)]dx= / f{x) exp[-G{x)]dx. 

From (4.2) we conclude that in each of the 13 cases in part (b) 

Pn,pn,K„{Sn/n^~^ e dx} =^ exp[-G(a;)] dx. 
This completes the proof of the theorem. ■ 

Two special cases of the scaling limits in Theorem 7.1 are worth pointing out. Given ^ > 
and A: 7^ 0, the sequence 

(3n = Pc and K,, = K{P,) - k/n' 

corresponds to the choice a = cxo in Theorem 7.1. With this sequence and with the same 
proofs, one obtains exactly the same limits as in cases 1, 3, 6, and 7 in this theorem with the 
same choices of 9, 7, and k. Similarly, given a > and b ^ 0, the sequence 

= log(4 - Vn") and Kr, = K{P,) 

corresponds to the choice ^ = 00 in Theorem 7.1. With this sequence and with the same proofs, 
one obtains exactly the same limits as in cases 1, 2, 4, and 5 in this theorem with the same 

choices of a, 7, and b. 

This completes our analysis of scaling limits for Sn/n^"^ with respect to Pn.i3„,Kn, where the 
sequence (/?„, Kn) converging to (/3c, Kc{/3c)) is defined in Theorem 7.1. In the next section we 
study MDPs for Sn/n}-''^ for appropriate sequences (/3„, K^) converging to (/5, K) e AliBVJC 
and for appropriate choices of 7. We obtain 1 MDP for (/3, K) e A,A MDPs for (/3, Kc{l3)) e 

B, and 13 MDPs for {p^, KdPc)) e C. 

8 18 MDPs for {/3n, Kn) ^ {/3, K) e AU B U C 

In this section we turn to a new problem, which is to formulate MDPs for Sn/n^^"' with respect 
to Pn,p„,K„, first for appropriate sequences {(3n, Kn) converging to {/3,Kc{(3)) G B, then for 
{Pn:Kn) Converging to {P,K) e A, and finally for (Pn^K^) converging to {Pc^K^iPc)) e 

C. These results are stated, respectively, in Theorem 8.1, Theorem 8.2, and Theorem 8.3. In 
proving the first result, we introduce the methods that are also used to prove the third. The 
proof of the MDP when (/3„, Kn) {P,K) e A proceeds differently from the proofs of 
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the other MDPs in this section, relying on the Gartner-Ellis Theorem. After the proof of that 
MDP, we will remark on why the same method cannot be used to prove all the MDPs in this 
section. Although an MDP is an LDP, we shall follow the example of [14], who in their study of 
Curie- Weiss -type models speak about an MDP whenever the exponential speed a„ of the large 
deviation probabilities satisfies a„/n ^ as n — > oo. Also see [12, §3.7]. 

When (Pri, Kn) Kc{(3)) e S U C, we will prove the MDPs by a method that is 

closely related to the proofs of the scaling limits earlier in this paper. Thus, rather than focus on 
the large deviation probabilities directly, we prove that Sn/n^^"' satisfies an equivalent Laplace 
principle. Despite the similarity in the proof of the scaling limits and the Laplace principles, the 
proof of the latter is much more delicate, requiring additional estimates not needed in the proof 
of the former. 

We start by considering the MDPs when (/?„, Kn) converges to (/?, Kc{(3)) E B. In order 
to formulate these limit theorems, we adapt the methods used in section 6, where we proved 
scaling limits for such sequences {/3n,Kn). For /3 G (0, /?c) let /5„ be an arbitrary positive 
sequence that converges to /3. Given ^ > and A; 7^ 0, we then define Kc{/3) as in (6.2). 

With this choice, part (b) of Theorem 4.3 implies that for any 7 > and R > and for all 
X e ]R satisfying Ix/ii^] < R there exists ^ e [—x/rf, x/rf] such that [see (6.3)] 

1 r^^^ 1 g'^'^^ (0) 1 

nG,^,KMn-) = + ,4.-1 + ^B^(a^M)x'. (8.1) 

The error terms Bn{^{x/n')) are uniformly bounded over n e N and x G {—Rn'^,Rn'), 

cP ^ 2/3, and <Vn(0) ^ ^f^M > 0. 
Given 7 e (0, 1), we define 

v = min{27 + ^- 1,47-1}. (8.2) 
In Theorem 6.1 we prove that when ^; = 0, Sn/'n}~'^ satisfies the scaling limit 

Pn,i3n,K„{Sn/n^~'^ G dx} =^ exp[- G {x)]dx , 

where 

G{x) = d{v, 2^ + 9- l)kpx^ + 6{v, 47 - l)c4X^ 

and C4 is defined in (6.5). As enumerated in Table 6.1, the 4 different forms of the limiting 
density depend on the values of 7 and 9 and the sign of k. 

In Theorem 8.1 we prove the analogous results on the level of MDPs. Assume that the 
quantity v defined in (8.2) is negative. Then, when {l3n,Kn) is chosen as in Theorem 6.1, 

Sn/n^~'^ satisfies the MDP with exponential speed n~'" and rate function r(a;) = G{x) — 
infj,gKG(y), where G is defined in the last display. We prove the MDP in Theorem 8.1 by 
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showing that when v < 0, Sn/n} satisfies the Laplace principle with speed n and rate 
function F; i.e., for any bounded, continuous function t/j 

1™ ^ log / exp[n"'' ^(5„/n^"^)] dPn,0„,Kn = sup{^(a;) - T{x)}. 

By Theorem 1.2.3 in [13] the fact that Sn/n}"'^ satisfies the Laplace principle implies that 
Sn/n^^'' satisfies the LDP with the same speed n~'" and the same rate function F; i.e., for any 
closed subset F in M 

limsup ^logP„,^„,/^„{5„/n^~^ e F} < - inf r(a;) 

and for any open subset $ in R 

liminf ^logP„,^„,;,J5„/n^-^ G $} > - inf F(a;). 

F is obviously a rate function. One easily checks that in all 4 cases given in part (b) of Theorem 
8.1 —V < 1. Hence n~'"/n — >^ as n — cxd, and so we have an MDP. In cases 1, 2, and 3, we 
have infj/giR G{y) = and thus F = G; in case 4, Inf^gK G{y) < 0. 

As in the scaling limits in Theorem 6.1, the rate function in Theorem 8.1 takes the 4 forms 
enumerated in cases 1, 2, 3, and 4 in Table 8.1. In case 2 the requirement that G{x) — > oo as 
|x| — > oo forces A; > 0. By contrast, in case 4, < is allowed. In case 1 we can also choose 
k to be any real number; this affects only the definition of the sequence Kn, not the form of the 
rate function. 

The forms of the rate functions reflect the influence, respectively, of B, of A, and of A and B. 
In each case the particular set or sets that influence the form of G depend on the speed at which 
{Pn: Kn) approaches Kc{P)) and the direction of approach. Case 2, which corresponds to 
the influence of A alone, has two subcases, labeled 2a and 2b in Table 8.1. 

In Figure 5 and in Table 8.1 we indicate the subsets of the positive quadrant of the O-j 
plane leading to the 4 cases of the MDPs in Theorem 8.1. Subcases 2a and 2b correspond, 
respectively, to the left half and the right half of the triangle labeled A in Figure 5. An interesting 
connection between the MDPs in Theorem 8. 1 and the scaling limits in Theorem 6. 1 is revealed 
by comparing Figure 5 with Figure 4, which exhibits the subsets of the positive quadrant of the 
O-'j plane leading to the 4 cases of the scaling limits in Theorem 6.1. The subsets labeled A, 
B, and A + B in Figure 4 are each a subset of the boundary of the set having the same label in 
Figure 5. The relevant boundaries in Figure 5 are labeled d'^A, d^B, and d'^{A + B), the first 
two of which are indicated by dotted lines. This relationship between the two figures is not a 
surprise because the sets labeled A, B, and A + B in Figure 4 are determined by solving v — 
while the sets having the same labels in Figure 5 are determined by solving v < 0. 
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Figure 5: Influence of B and on MDPs when (/3„, K„) (/?, G B 



Theorem 8.1. For fixed /3 G (0, Pc), let (3n be an arbitrary positive sequence that converges to 
(3. Given 9 > Q and k 0, define 



= K{(3n) - k/n\ 

where K{l3) = (e^ + 2)/{4(3)for /? > 0. Then Kn) K,{(3)) G B. Given 7 G (0, 1), 



The following conclusions hold. 

(a) Assume that v = min{27 + 6* — 1, 47 — 1} satisfies v < 0. Then with respect to Pn,i3„,Kn' 
SnI n^'^ satisfies the Laplace principle, and thus the MDP, with exponential speed n~'" and rate 
function T{x) = G{x) — infj^gK G{y). 

(b) We have v < if and only if one of the 4 cases enumerated in Table 8.1 holds. Each 
of the 4 cases corresponds to a set of values of'y and 9, a choice of sign ofk, the infiuence of 
one or more sets B and A, and a particular exponential speed and a particular form of the rate 
function in part (a). The function G appearing in the definition of the rate function is shown in 
column 5 in Table 8.1; in case 4 the nonzero constant infj^g^ in the definition of the rate 
function is not shown. In case 1 the choice ofkEM. does not affect the form of the rate function. 



we also define 



G{x) = 6{v, 2-f + 9- l)kl3x^ + 6{v, 47 - l)cix\ 



(8.3) 



where C4 > is given by 




2[2(3K,{(3)]\4-e'') _ (e^ + 2)^(4 - e^) 
4!(e/5 + 2)2 ~ 23-4! 
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case 


values of 7 


values of 


exp'l 


function G in 


influence 






speed 


rate function F 


1 


7e(0,i) 


e>2-f 


^1-47 


4 


B 








C4 > 0, A; e M 


2a 


7 e (0, i] 


(0,27) 


^1-27-^^ 




A 








A; > 


2b 




^ e (0, 1 - 27) 


^l-27-y 




A 








k>0 


3-4 


7e(0,^) 


e = 27 


^1-47 


k(3x'^ + CiX^ 


A + B 








k^O 



Table 8.1; Values of 7 and 9, exponential speeds, and rate functions in part (b) of Theorem 8.1 

Proof. We first prove part (b) from part (a) and then prove part (a). 

(b) We have v < in the following 4 mutually exclusive and exhaustive cases. As (8.3) 
makes clear, f = 47 — 1 < corresponds to the influence of B and v = 27 + 6* — 1 to the 
influence of A. 

• Case 1: Influence of B alone, -u = 47 - 1 < 0, 47 - 1 < 27 + ^ - 1, and A; G M. In 
this case 7 G (0, 1/4) and 9 > 27, which corresponds to the second and third columns 
for case 1 in Table 8.1. 

• Case 2: Influence of A alone. ?; = 27 + 6* - 1 < 0, 27 + 6' - 1 < 47 - 1, and A; > 0. 
In this case < 9 < min{27, 1 - 27}. Since < 27 < 1 - 27 7 G (0, 1/4] and 
0< 1 — 27<27-^7G (1/4, 1/2), case 2 corresponds to the second and third columns 
for case 2a and case 2b in Table 8.1. 

• Cases 3-4: Influence of A and fi. f = 47 — 1 = 27 + 6^ — 1 < 0, A; > for case 3, and 
A; < for case 4. In these cases < 7 < 1/4 and 9 — 27. Hence case 3-4 correspond to 
the second and third columns for cases 3-4 in Table 8.1. 

In cases 1, 2, 3, and 4 we have, respectively, G{x) = c^x"^, G{x) = kjSx^ with A; > 0, G{x) = 
kPx"^ + c^x^ with A; > 0, and G{x) = k(3x'^ + c^x^ with A; < 0. In combination with part (a), 
we obtain the 4 rate functions given in the last column of Table 8.1. 

(a) Our strategy is to prove that with respect to Pn,f3„,K„ x Q, Sn/n^~'^ + Wn/n^^"^'^ sat- 
isfies the Laplace principle with exponential speed n~'" and rate function V. In order to prove 
the Laplace principle for Sn/n}~^ alone, we need the following estimate, which shows that 
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Wn/n}^'^ is superexponentially small relative to exp(n for any S > 

limsup^logQ{|W^„/n^/2-^| > S} ^ -oo. (8.4) 

n— »oo ^ 

According to Theorem 1.3.3 in [13], if with respect to Pn,i3n,K„ x Q, Wn/n^^"^^"' + Sn/n^~"' 
satisfies the Laplace principle with speed n^'" and rate function T, then with respect to Pn,(3n,K„, 
Sn/n}~'^ satisfies the Laplace principle with speed n~'" and rate function F. Since the Laplace 
principle implies the MDP [13, Thm. 1.2.3], part (a) of the present theorem will be proved. 

We now prove (8.4). Denote the variance (2/3„A'„)^^ of Wn by af^. Since and Kn are 
bounded and uniformly positive over n, the sequence cr^ is bounded and uniformly positive over 
n. We have the inequality 

Q{\WJn"''-'\ > S} = Q{\N{0,al)\>n"^-'S} 

Hence (8.4) follows if 1 — 27 > —v. Since 7 and 9 are both positive, this is easily verified to 
hold when either f = 47 — 1 or = 27 + ^ — L 

We now turn to the Laplace principle for Sn/n^^"' + Wn/n^^'^^'^ ■ Let '(jj be an arbitrary 
bounded, continuous function. Choosing / = exp[n^^t/'] in Lemma 4.1 yields 



/ exp 



^1-7 ^1/2-7 
1 



d{Pn,p„,K^ X Q) (8.5) 
expln'^'i/jix) - nGf3^,Kr,{x/n^)] dx, 



L^M-nG/3„,K„{x/n'r)] dx 

In order to obtain the appropriate expansion of nG p^^x^x / n^) in this display, we multiply the 
numerator and denominator of the right hand side of (8.1) by n~^, obtaining 

where 

r \ '* 2 , 1 ^pl,Kn(^) 4 , 1 c ft-f / "/w 5 
Gn(X) — „ ,■■ —X H — -. ^^-v; X H — = — - — BniEix rfux . 

The proof of the Laplace principle for Sn/'n}~'^ + Wn/n}^'^''^ rests on the following prop- 
erties of nGp^^Knix/n"') = n~'"Gn{x), which in turn are consequences of the Taylor expansion 
of Gn{x) just given. Because of the estimate (8.4) on Wn/'n}-/'^~'^ , the inequality in (8.6), and 
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the uniform convergence of Gn to G expressed in item 3 below, the proof of the MDPs, though 

analogous, is more delicate than the proof of the scaling limits in section 6, for which the a.s. 
convergence of Wn/n^^'^~'^ to 0, the pointwise convergence of G p^^K„{x / ni) to G{x), and the 
lower bound (6.7) suffice. 

1. There exists /? > and a polynomial H with the properties that H{x) ^ oo as — > oo 
and for all sufficiently large n and all x G M satisfying \x/rO\ < R 



In case 1 when A; > as well as in cases 2 and 3, H{x) = G{x)/2; in case 1 when A; < 
and in case 4, which corresponds tok < 0, H{x) — —2\k\Px'^ + 04x^/2. 

2. Let A = sup^^^{ip{x) — G{x)}. Since H{x) — > 00 and G{x) — > 00, there exists M > 
with the properties that 



the supremum of ^ — G on R is attained on the interval [— M, M], and the supremum of 
— G on M is attained on the interval [— M, M]. In combination with item 1, we see that 
for all n e N satisfying Rn'^ > M 



3. Let M be the number selected in item 2. Then for all x e R satisfying \x\ < M, 
Gn{x) — 'n}-^'"Gf}^^Kn{^/'^'^) converges uniformly to G{x) as n — > 00. 

Since nG p^^Xr^i^ / n'^) — n'''"Gn{x), item 3 implies that for any 5 > and all sufficiently 
large n 



nGp„,K„{x/n^)>n-''H{x). 



sup {ij{x) - H{x)} < -\A\ - 1, 



|a;|>M 



sup {n'^^ix) - nGp^,KMn'^)} < -^"'(1 A| + 1). 



(8.6) 



M<\x\<Rn'< 




In addition, item 2 implies that 




exp[n-'"ij{x) - nGp^,Kr.{x/n^)] dx < 2i?n^ exp[-n"^(| A| + 1)]. 
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Since ip is bounded, the last two displays show that there exist > and ae e R such that for 
all sufficiently large n 

J{|a;|<RnT} 

Since —v G (0, 1), we conclude from part (d) of Lemma 4.4 the existence of ay > such that 
for all sufficiently large n 

{\x\>Rn^} 

We now put these three estimates together. For all sufficiently large n we have 

exp{—n~^5) / ex.'p[n~'"{ilj{x) — G{x))]dx 

J{\x\<M} 

< / exp[n~''ip{x) - nGi3^^K„{x/n'^)]dx 



< exp(n-^5) / exp[n-''(^(a;) - G{x))] dx + 5^, 



'{\x\<M} 

where 

Sn < 2i?n^exp[-n-^(|A| + 1)] +2a5exp(-r^a7 + r^-^||V'||oo)■ 
Since —v<l and since by item 2 



lim ^^log / exp[n '"{'i/j{x) — G{x))]dx 

n-voon ^ J{\x\<M} 

— sup {ipix) — G{x)} = sup{'^(a;) — G{x)}, 

\x\<M xeM. 



we have 



sup{ip{x) — G{x)} — S 



< 



< 



liminf^log [ exp[n ''^(x) - nG'^„,;^„(x/n^)] 
lim sup — log / exp[n-'''ijj{x) - nG p^^xS^ / n^)] dx 



< sup{ilj{x) - G{x)} + S, 
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and because 5 > is arbitrary, it follows that 

lim log / exp[n~'"ip{x) — nGa dx — sup{ip{x) — G{x)}. 

Combining this limit with the same limit for ^ = 0, we conclude from (8.5) that 



lim log 



/ exp 









s„. ^ w„. 



^1-7 nV2-7 



sup{ip{x) — G{x)} + inf G{y) — sup{ip{x) — r(x)}. 



This completes the proof that with respect to Pn.8n,Kn x Q, Sn/n^^^ + Wn/n^^"^^^ satisfies 
the Laplace principle with exponential speed and rate function V. Since Wn/n^^'^~'^ is 
superexponentially small, we obtain the desired Laplace principle for Sn/n^~'^ with respect to 
Pn,i3n,Kn- The proof of the theorem is complete. ■ 

We next formulate the MDP for Sn/n^^'^ when (/5„, Kn) is an arbitrary positive sequence 
that converges to K) e A; thus /? and K satisfy < /3 < /3c and < < Kc{l3). Because 
in this case the normal random variable Wn contributes to the limit, we are not able to prove 
the MDP as we proved Theorem 8.1. Instead we use the Gartner-Ellis Theorem. The following 
theorem is also valid for (3 > jSc and < i^' < Kc{(3), and the proof is essentially the same. 
The key observation is that for p > (3^, we have K{(3) = {e^ + 2)/(4/5) > Kc{P) [22, Thm. 
3.8]. Hence if K < K^{(5), then also K < K{(3) and thus ^^^(O) in (8.7) is positive. 

Theorem 8.2. Let Kn) he an arbitrary positive sequence that converges to {P, K) e A. 
Let 7 be any number in (0, 1 /2). Then with respect to Pn,f3„,Kn' ^n/ 'n^^'"' satisfies the MDP with 
exponential speed n^^^^ and rate function j3[K[j3) — K]x^. Thus the limit is independent of the 
particular sequence {f3n, Kn) that is chosen. 

Proof. For n e N and t G M we use the monotone convergence theorem to replace / in Lemma 
4.1 by exp(n^~'^'^tx). We then use the Taylor expansion in part (a) of Theorem 4.3 and the fact 
that G^^^^i(^{0) given in (4.9) converges to 



(,) _ 2/3X[X(/3) - K] 
G^M^) , (8.7) 



which is positive since < K < Kc{P) — K{/3). As in the proof of part (a) of Theorem 8.1, 
there exists M > such that the supremum of tx — G^^^(0)a:^/2 is attained on the interval 
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[— M, M] and the following calculation is valid 

J An 



n—*oo n 



1-27^ / 



n 



1-7 



+ 



n 



1/2-7 



lim log / exp[n^ '^"'tx - nGp^^K„{x/ri')] dx 

n— >oo n ' Jjj 



— lim 



1 



n-+oo 77, 



1-27 



log / exp[-nG';3„,;^„(a:;/n^)]dx 



= lim 

n— >oo 7T, 



^ log / exp \n'-'^ (tx - Gf^{0)xy2) 
lim log / exp [-Ggi,(0)xV2] 



sup {te - 4'5^(0)a;V2| + inf |G'^'i^(0)a;V2| 



{|x|<M} 



2^0) ' 



Since VF„ is an iV(0, (2/3„-ft'„) ^) random variable and is independent of S„ 



lim 

n— >oo 77, 



1-27 

lim 

n— >oo 77, 



log / exp 



A" 



71 



-27^. 



exp 



771-T 
1 



A"xn 



71 



n 



1-7 



1/2-7 



71 



• lim ^Vlog / exM/^~^tWn\ dQ 



The Gartner-Ellis Theorem [16] now implies that Sn/n}"'^ satisfies the MDP with exponential 
speed 'n}~'^'^ and rate function 



)< — — - • 



/(x)-sup^.. 2 2/3TO)-X] 



This completes the proof. 



In the context of the proof of the preceding theorem, it is worthwhile pointing out that the 
Gartner-Ellis Theorem cannot be used to prove all the other MDPs for Sn/'n}~'^ in this section. 
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For example, consider the MDPs in Theorem 8.1. For any t eM. one calculates 



g{t) = lim ^log / exp 



n t\ — r— + 



^1-7 72,1/2-7 



= sup{tx - G{x)} + inf G{y) = sup{tx - [G{x) - G]}, 



where (5 = infj^gR G{y). Thus equals the Legendre-Fenchel transform of G — G. If G — G 
is strictly convex on M, as it is in cases 1, 2, and 3 in Theorem 8.1, then g is differentiable on M 
[32, p. 253]. Hence by the Gartner-Ellis Theorem, Sn/n^~"' satisfies the MDP with exponential 
speed n^'" and rate function given by the Legendre-Fenchel transform of g, which is G — G. 
In cases 1, 2, and 3 in Theorem 8.1, (5 equals 0, and we recover the form of the rate function 
in column 4 of Table 8.1. However, the situation is different in the MDP in case 4, in which 
G{x) = kSx^ + c^x'^ with /c < 0. Here G < 0, G is not convex on all of M, and g is not 
differentiable on M. As a result, the Gartner-Ellis Theorem cannot be applied to obtain the lower 
large deviation bound for all open sets and thus to obtain the MDP. In addition, the Legendre- 
Fenchel transform of g equals on a symmetric interval containing the origin, and thus it does 
not coincide with G G on this interval. A similar situation holds in Theorem 8.3, in which 
we derive 13 MDPs for suitable sequences Kn) ~^ (Pc, Kc{Pc))- In cases 1-4, 6, 8, and 10 
in that theorem, the coefficients in the polynomial G are all positive, and so G is strictly convex 
and G = 0. Hence the corresponding MDPs can be derived via the Gartner-Ellis Theorem. 
However, in all the other cases except for case 12 with k sufficiently large, the polynomial G 
is not convex on all of M; as in case 4 in Theorem 8.1, the Gartner-Ellis Theorem cannot be 
applied to obtain the MDP. 

We now consider the final class of MDPs in this section. This class arises when Kn) 
converges to (/3c, i^c(/^c)) along the same sequences considered in Theorem 7.1, where we 
proved scaling limits for Sn/n^~'^ for 7 e (0, 1/2). Given a > 0, 9 > 0, b 0, aad k 0, 
these sequences are defined by 

(3n = log(4 - Vn") = log(e'5= - 6/^") and = K{[5n) - k/n'. (8.8) 

For these sequences the parameter that plays the role of v in Theorem 8.1 is 

w — min{27 -|- ^ — 1, 47 -|- a — 1, 67 — 1}. 

The 13 forms of the scaling limits of Sn/'n}~'^ are proved in Theorem 7.1 under the assumption 
that w — 0. We now assume that w < 0. Using the same Taylor expansion that was used to 
deduce these scaling limits [Thm. 4.3(c)], one deduces the 13 forms of the Laplace principles 
for Sn/n^~'^ . These Laplace principles and the equivalent MDPs are stated in the next theorem 
along with the choices of 7, a, b, 6, and k leading to the 13 forms of the rate function. The only 
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requirement on b and k is that G{x) — > oo as |x| — > oo. This requirement forces ft > in case 
2 and A; > in case 3. The proof of the MDPs in the next theorem is omitted because it follows 
the same pattern of proof of Theorem 8.1. 

As in Theorem 7.1, there are further possibilities concerning the sign of b and k. In all the 
cases in which no term appears in the scaling limit (cases 1, 3, 6, 7), we can choose either b 
to be any real number. Similarly, in all the cases in which no x"^ term appears in the scaling limit 
(cases 1, 2, 4, 5), we can choose either k to be any real number. Although the choice of 6 or /c 
affects the definition of the sequence Kn), it does not affect the form of the rate function. 

Theorem 8.3. Given a > 0, 9 > 0, b 0, and k 0, consider the sequence (/3„, Kn) defined 
in (8.8). Then Kn) {Pc, Kc{Pc))- Given 7 e (0, 1), we also define 

G{x) = S{w, 2j + 9- l)kpcx^ + S{w, 47 + a - l)bc^x^ + S{w, 67 - l)c6X^, 

where C4 = 3/16 and cq = 9/40. The following conclusions hold. 

(a) Assume that w = min{27 + ^ — 1, 47 + a — 1, 67 — 1} satisfies w < 0. Then with respect 
to Pn,)3n,Kn' Sn/n^~'^ Satisfies the Laplace principle, and thus the MDP, with exponential speed 
n~'^ and rate function T{x) — G{x) — inf j^gu G{y). 

(b) We have w < Q if and only if one of the 13 cases enumerated in Table 8.2 holds. Each 
of the 13 cases corresponds to a set of values of'y, a, and 9; a choice of signs ofb and k; the 
influence of one or more sets C, B, A; and a particular exponential speed and a particular form 
of the rate function in part (a). The function G appearing in the definition of the rate function is 
shown in column 5 in Table 8.2; when inf^gK G{y) ^ 0, this additive constant in the definition 
of the rate function is not shown. The form of the rate function is not affected by the choice of 
b ^Win cases 1, 3, 6, and 1 and by the choice ofk&R. in cases 1, 2, 4, and 5. 
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case 


values of 7 


values of a 


exp'l 


function G in 


influence 




values of d 


speed 


rate function F 


1 


7e(0,i) 


a > 27 


^1-67 




C 




^ > 47 




C6 > 0, 6 e M, A; e M 


2 


7e(0,i) 


a e (0,min{27,l -47}) 


^1-47-a 




B 




6* > 27 + a 




6 > 0, C4 > 0, e R 


3 


7e(0,i) 


e e (0,min{47,l - 27}) 






A 




a > max(^ - 27, 0) 




k>0,beR 


4-5 


7e(0,^) 


a = 27 


^1-67 




B + C 




^ > 47 




6 ^ 0, A; e M 


6-7 


7e(0,i) 


a > 27 


^1-67 




A + C 




^ = 47 




A; ^ 0, 6 G K 


8-9 


7=(0,i) 


a e (0,min{27, 1 - 47}) 


^1-47-a 




A + B 




6* = 27 + a 




A; ^ 0, 6 > 


10-13 


7e(0,i) 


a = 27 


^1-67 




A + B + C 




^ = 47 




A; ^ 0, 6 7^ 



Table 8.2; Values of 7, a, and ^, exponential speeds, and rate functions in part (b) of Theorem 8.3 



As discussed in section 2, the MDPs listed in Table 8.2 yield a new class of distribution 
limits for Sn/n^~'^ in those cases in which the set of global minimum points of G contains 
nonzero points. These are the cases in which the coefficients of G are not all positive: cases 
5{h < 0), 7 (A; < 0), 9 (A; < 0), 11 (A; < 0, 6 > 0), 12 (A; > 0, 6 < 0), and 13 {k < 0, 
h < 0). In all these cases except for case 12, we obtain the limit (2.13). Case 12 exhibits the 
most complicated behavior, giving rise to the limit (2.14) for the critical value k = 56^/[2''/5c]. 
These limits and the underlying physical phenomena are now being investigated for a class of 
non-mean-field models, including the Blume-Emery-Griffiths model [19]. 

This completes our study of limit theorems for the BEG model in the neighborhood of the 
tricritical point {(3^, K^{Pc)) E C, in the neighborhood of second-order points (0, K,{(3)) G B, 
and in the neighborhood of single-phase points K) G A. It is an unexpectedly rich and 
fruitful area of research, one that we hope will inspire similar investigations for other statistical 
mechanical models. 
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